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unsteady  flow.  Instead,  they  independently  proposed  later  the  so- 
called  MRS  separation  criterion  (the  velocity  cross  section  near 
the  corresponding  separation  point  is  called  the  MRS  cross  section). 
Then,  Sears  and  Telionis(197l) (4)  pointed  out  that  solution  to  the 
unsteady  boundary-layer  equation  exhibited  a  Goldstein- type  sin¬ 
gularity  at  separation.  From  that,  the  "singular  point  separation 
criterion"  was  proposed.  Vidal( 1959) (5)  and  Ludwig( 1964) (6) 
verified  the  MRS  criterion  in  a  steady  flow  experiment  on  a  moving 
wall.  The  MRS  criterion  and  the  singular  point  separation  criterion 
were  also  verified  under  the  downstream-moving  wall  conditions  (or 
the  wall  did  not  move,  but  the  separation  point  moved  upstream)  by 
Telionis  and  Werle( 1973) (7)  in  numerical  solutions  to  the  steady 
boundary-layer  on  a  downstream-moving  wall;  by  Telionis,  Tsahalis 
and  Werle( 1973) (8)  with  the  difference  solution  to  the  unsteady 
Howarth's  linear  deceleration  boundary-layer;  and  by  Williams  and 
Johnson( 1974) (9) ,  Williams( 1981 ) ( 10) ,  and  C.  T.  Wang  and  S.  F.  Shen 
(1978) (11)  in  the  semisimilar  solution  to  the  Howarth's  linear 
deceleration- type  and  Falkner-Skan- type  unsteady  boundary-layers. 

In  the  case  of  upstream-moving  walls,  however,  the  MRS  criterion 
was  not  verified  by  Danberg  and  Fansler(l2)  with  the  similar  solu¬ 
tion  to  the  steady  boundary  layers;  by  Tsahalis(l3)  in  the  numerical 
solution  to  the  unsteady  boundary  layers;  and  by  Inoue(l4)  with  the 
numerical  solution  to  the  simplified  NS  equation.  Recently,  Dommele-n 
and  Shen  calculated  the  instantaneously-starting  boundary-layer  flows 
on  a  cylindrical  surface  using  Lagrange's  method.  They  found  that 
when  the  cylinder  moved  from  the  starting  point  through  a  distance, 
0.75  of  its  diameter,  a  singular  point  appeared  at  110°  from  the 


front  stagnation  point,  then  moved  forward  along  the  surface  at  a 
speed  0.52  times  that  of  the  cylinder,  and  finally  stopped 
at  104.5°  steady  position.  They  also  found  that  although  the 
criterion  could  serve  at  the  singular  point,  the  friction  stress 
curves  were  smooth  and  even,  while  the  displacement  thickness 
increased  dramatically.  Therefore,  the  singular  point  was  not  of 
the  Goldstein- type ,  which  was  controlled  by  the  viscous  friction 
force  (also  called  a  viscous  singular  point),  but  of  another  type, 
which  was  called  inviscid  singular  point.  K.  C.  Wang's  unsteady 
similar  calculation  has  confirmed  Dommelen  and  Shen's  results  to 
some  extent.  However,  according  to  Cebeci(l8),  Dwyer  and  Sherman 
(19),  within  a  limited  period  of  time,  there  is  no  singular  point 
in  an  instantaneously-starting  cylinder. 

In  1950's,  research  work  was  conducted  on  the  origin  of 
separation  and  the  boundary  line  (separation  line)  patterns  in  the 
separation  region  for  three-dimensional  steady  flows.  To  date,  a 
lot  of  theories,  experimental  work(20-27),  and  review  papers( 28-30) 
are  available  on  the  topic.  However,  Brown. and  S tewartson( 28) ,  and 
Wiiliams(30) ,  et  al.  pointed  out  that  there  still  existed  a  lot  of 
controversy  on  this  subject.  For  example,  about  the  characteristics 
of  the  separation  line,  Lichelbrenner  and  Cudart(31);  Faskell ( 20) ; 
and  K.  G.  Wang  believed  ealier  that  the  separation  line  was  an  en¬ 
velope  of  the  limiting  streamlines  on  a  wall.  However,  according  to 
Legendre( 21 , 32) ,  Lighthill ( 22)  (who  developed  the  concept  of  fric¬ 
tion  force  lines  which  are  usually  considered  to  be  equivalent  to 
the  limiting  streamlines)  and  Hunt(33);  and  Tobak  and  Feak e(34-*0, 
a  separation  line  is  a  limiting  streamline  on  the  wall  and  excent 


for  the  singular  point,  it  is  impossible  for  the  nearby  limiting 
streamlines  to  intersect  the  separation  line  which  is  a  "convergent 
asymptote"  within  a  limited  distance.  Also,  Brown  and  Stewartson( 28) 
questioned  the  possibility  for  the  separation  line  of  the  flow  des¬ 
cribed  in  the  boundary-layer  equation  to  exhibit  a  Gold  stein- type 
singularity  similar  to  the  separation  point  of  the  two-dimensional 
unsteady  flow.  Recently,  Cousteix’s  work  using  approximate  integral 
relation  techniques  shows  the  possibility  of  exhibiting  such  a  sin¬ 
gularity. 

The  origin  and  terminal  of  the  three-dimensional  seoaretion 
line  are  other  interesting  topics.  Legendre  (21)  and  Lighthill  (22) 
suggested  that  the  origin  of  the  separation  line  should  be  a  saddle 
point,  while  Kaskell(20),  K.  G.  Wang(24),  et  al.  found  through  their 
experimentation  and  calculations  that  there  were  two  types  of  sepa- 
ation  origin:  bubble  and  vortex-layer  separation,  as  defined  by 
iXaskell;  open-  and  closed- type  separation,  as  defined  by  K,  C.  Wang; 
and  the  integral  and  partial  separation,  as  recently  defined  by 
Tobak,  Peake,  et  al.  For  bubble  or  closed- type  (or  integral)  seoa- 
ration,  the  origin  of  the  separation  line  is  actually  a  saddle  sin¬ 
gular  point,  but  the  origin  of  vortex  or  open-type  (or  partial)  sep¬ 
aration  is  not  a  singular  point.  Research  on  the  termination  of  the 
separation  line  shows  that,  exceDt  for  the  singular  point  of  zero  fric¬ 
tion  stress,  the  separation  line  chancres  continuously,  and  when  it  meets 
a  focus  or  a  node,  the  line  terminates  at  such  a  singular  point. 
Because  of  the  connection  of  the  singular  point  with  separation, 
legendre( 21 ) ,  Oswatitsch( 23) ,  and  Lighthill( 22)  studied  the  singular 
point  distribution  on  a  material  surface  using  topological  techniques. 


Hunt(33),  et  al.  extended  the  material  surface  singular  point 
and  topological  analysis  to  include  body  cross-section  in  the  out¬ 
side  flow  field,  and  proposed  semi-nodal  and  semi-saddle  ooint  con¬ 
cepts.  Tobak  and  Peake(34)  further  investigated  the  applicability 
of  topology  to  the  three-dimensional  separation  flow  and  put  for¬ 
ward  topological  structure  and  structure  stability  theories( 35 , 36) . 
ihe  authors(38)  and  Liu  Koujie,  et  al .  also  anslyzed  these  singular 
points. 

In  summary,  the  following  questions  concerning  the  origin  of 
separation  and  the  patterns  of  the  two-dimensional  unsteady  and 
three-dimensional  steady  flows  need  to  be  answered. 

(1)  For  two-dimensional  unsteady  flow,  are  KRS  criterion  and 
the  singular  point  separation  criterion  universal  and  equivalent  to 
each  other?  Does  the  singular  point  given  by  the  singular  point 
separation  criterion  exhibit  viscous  or  inviscid  behavior?  Does 
the  separation  singular  point  appear  after  a  limited  duration  of 
time? 

(2)  For  three-dimensional  steady  flow,  what  are  the  seoaration 
conditions?  Is  the  separation  line  an  envelope  of  the  wall  limiting 
streamlines  or  an  asymptote  of  the  limiting  streamlines  themselves 
and  nearby  ones?  Does  the  separation  line  for  the  flow  described  in 
the  steady  boundary  -layer  equation  exhibit  a  singularity?  How  does 
the  separation  line  originate  and  terminate?  What  are  the  singular 
point  distribution  and  the  separation  line  patterns  at  the  singular 
point? 

Koreover,  according  to  K.C.  Wang's  analogue  seoaration  criter¬ 
ion,  under  certain  conditions,  the  two-dimensional  unsteady  flow  can 


be  analogized  as  the  three-dimensional  steady  flow,  and  the  seoa- 
ration  criteria  for  the  two-dimensional  unsteady  flow  and  three- 
dimensional  steady  flow  should  be  transformed  from  each  other. 

It  is  worth  investigating  if  a  transformation  of  the  separation 
criteria  now  available  for  the  two  types  of  flows  can  be  performed. 

The  objective  of  this  paper  is  to  discuss  and  review  the 
above-mentioned  problems,  and  to  introduce  the  author's  work  on 


the  topic. 


2.  TWO-DI EENSIONAL  UNSTEADY  FLOWS 


Sears,  koore  and  Kartunian  pointed  out  the  similarities  be¬ 
tween  unsteady  flow  over  fixed  walls  and  steady  flow  over  moving 
walls.  Sears  defined  separation  of  two-dimensional  unsteady  flow  as 
"the  appearance  of  the  stagnation  point  and  the  separation  line 
which  separates  the  boundary  layer  and  the  wake  in  a  coordinate  sys¬ 
tem  moving  with  the  separation  point"(see  literature  (3-13)). 
According  to  this  definition,  a  separation  point  is  of  stagnation 
type.  Toore,  Rott  and  Sears  have  also  provided  the  flow  patterns 
near  the  separation  point  (RRS  sketches),  as  shown  in  Fig.  1.  In 
the  figure,  (a)  corresponds  to  steady  flow  over  fixed  walls,  and 
(b)  and  (c)  correspond  to  the  moving  wall  condition.  After  a  care¬ 
ful  analysis,  O'Brien  pointed  out  in  1981  that  unsteady  laminar 
separated  flow  exhibited  three  patterns  which  were  similar  to  these 
of  a  steady  flow,  as  shown  in  Fig.  2.  When  only  the  vicinity  of 
the  stagnation  point  is  considered,  three  O'Brien's  patterns  are 
completely  consistent  with  the  MtS  sketches.  lost  studies  on  two- 
dimensional  unsteady  separated  flows  today  start  with  this  defi¬ 
nition  of  separation  and  these  separation  sketches. 


m 
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In  the  two-dimensional  unsteady  separated  flow  studies,  one 
of  the  important  topics  is  criteria  for  separation.  The  following 
are  those  suggested  in  the  literature. 

(1)  MRS  (Moore,  Rott  and  Sears)  criterion( 1-3) .  Conditions 


for  separation  to  appear  are  that  at  the  separation  noint 
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where  u  is  the  velocity  component  in  x  direction;  y,  the  direction 
perpendicular  to  the  material  surface;  and  xoy,  the  coordinate  sys¬ 
tem  moving  with  the  separation  point. 

(2)  O'Brien's  criterion(42) .  The  separation  conditions  are 


that  at  the  separation  point 
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(2.2) 


where  v  is  the  velocity  component  in  y  direction. 


(3)  Singular  point  separation  criterion  (Sears  and  Telionis) 

(4).  The  separation  conditions  for  the  flow  described  in  the 


boundary-layer  equation  are 


0 


(2.3) 


where  x  is  the  partial  shear  stress,  namely,  the  separation  noint 
is  a  Goldstein- type  singular  point. 

(4)  Singular- type  parameter  criterion  (Fansler  and  Danberp) 
(43).  Condition  for  the  boundary-layer  separation  is  that  at  the 


senaration  ooint 


(&l~ 


(2.4) 


where  H  is  the  ratio  of  displacement  thickness  of  the  boundary  layer 
to  the  momentum  thickness;  and  k,  the  ratio  of  the  energy  thickness 
to  the  momentum  thickness. 

(5)  Analogue  separation  criterion  (K.  C.  Wang) (17)  states  that 


for  two-dimensional  unsteady  flow,  the  separation  criterion  can  be 
established  based  on  that  for  three-dimensional  flow. 


In  addition,  C.  T.  Wang  and  S.  F.  Shen(ll)  also  proposed  a 
separation  criterion  which  is  similar  to  MRS  criterion. 

Many  investigators  have  been  devoted  to  examining  and  evalu¬ 
ating  the  above  separation  criteria.  It  is  generally  accepted  that 
MRS  criterion  and  zero  vortex  criterion  are  valid  under  conditions 
with  flow  over  downstream-moving  walls,  but  lead  to  erroneous  re¬ 
sults  in  cases  of  flow  over  upstream-moving  walls;  zero  vortex 
criterion  is  not  valid  especially  when  the  streamlines  vassing 
through  the  separation  point  are  not  normal  to  each  other(42). 

The  singular  point  separation  criterion  and  singular- type  parameter 
criterion  can  only  be  applied  to  certain  boundary  layer  flows.  The 
analogue  separation  criterion  is  also  dependent  on  the  criterion 
for  separation  of  three-dimensional  steady  flow,  which  is  still 
under  discussion.  In  view  of  the  above-mentioned  facts,  we  studied 
the  separation  of  incompressible  viscous  flow  using  singular  point 
theory  for  differential  equations  based  on  Sears'  definition  of  sep 
aration,  and  give  the  following  conditions  for  separation( 40 , 4 1 ) : 

(1)  In  a  viscous  flow  field,  if  uo=vo=0  at  point  "n",  and 

'-feX'fcUwl  >•  (2-5) 


is  satisfied,  then  point  ’’o"  is  a  separation  point.  If  F<0,  ooint 
"o"  is  not  a  separation  point.  Here,  xy  is  the  coordinate  system 
moving  with  ”o",  and  u  and  v  are  velocity  comoonents. 

(2)  If  uo=vo=0  and  F=0  at  point  "o",  and  fi0=0  and 

(g1 +  #')/<?*  =  0  ( 
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and  At-  (**»/***)„  B,-^**/***),,  C,-(dV***y) «.  A“(*V»y*)» 


,  then  "o' 


is  a  separation  ooint. 


Based  on  this  study,  we  can  draw  the  following  conclusions 


(41): 


(1)  From  condition  (l),  there  exists  a  zero  u  line  (u=0)  mass¬ 


ing  through  the  separation  point.  This  is  corresponding  to  the  Mu- 


sketch  in  Fig.  1(b).  From  condition  (2),  there  are  two  zero  u  lines 


passing  through  the  separation  point,  which  correspond  to  the  FRJ 


sketches  in  Fig.  1(a)  and  1(c). 

(2)  If  the  HRS  conditions  hold  at  point  "o",  and  =?  0 

then  condition  (l)  must  be  satisfied  and  point  "o"  is  a  separation 


point.  Thus,  condition  (l)  contains  HUS  criterion.  In  this  case, 


one  of  the  slopes  for  streamlines  passing  through  "o"  is  infinity 


(here  (*V***y)«  is  not  an  infinity).  Another  slope  is-i-^-^-^  /("^r)« 

(3)  If  zero  vortex  criterion  holds  at  point  "o"  and 
or  JifO,  i  then  condition  (l)  must  be  satisfied  and  "o  "  must 

be  a  separation  point.  Thus  condition  (1)  contains  zero  vortex 


criterion. 


(4)  For  the  separated  flow  which  satisfies  condition  ( 1 ) , 


two  streamlines  passing  through  separation  point  "o"  are  not  normal 

/  av  \  (  =  0 

to  each  other,  except  that  dx  )t~\  ay  ’  . 

(5)  For  steady  flow  over  a  fixed  wall,  if  (-T(r)t,  “  0  ^ 0  ’ 


is  satisfied  at  point  !'o"  on  the  wall,  then  condition  (2)  is  satis¬ 


fied.  Thus,  point  "o"  is  either  a  separation  point  or  an  attach¬ 


ment  point.  This  is  the  Frandtl’ s  criterion  for  steady  flows. 


(C)  For  flow  descrioed  in  the  boundary  layer  equations,  if 


at  point  "o" 
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and  (*■/**),«(*?/**),*■  0 ,  ,  then  "o"  is  generally  a  Sold  stein- tyne  sin¬ 

gular  point,  and  condition  (2)  is  satisfied.  Therefore,  "o"  is  a 
separation  point,  and  the  slope  of  one  streamline  passing  through 
"o"  is  infinity.  This  is  singular  point  separation  criterion. 

(7)  For  the  separated  flow  satisfying  condition  ( 1 ) ,  MRS  and 
zero  vortex  conditions  are  not  necessary.  But,  it  is  necessary  for 
the  separated  flow  satisfying  condition  (2). 

The  above  conclusions  tell  us  the  following: 

(1)  The  separation  criterion  conditions  based  on  Sears' 
definition  divide  flows  into  two  categories:  one  with  a  zero  u  line 
and  another  with  two  zero  u  lines.  These  are  consistent  with  both 
KRS  and  O'Brien  sketches. 

(2)  The  PUS,  zero  vortex,  singular  point  separation  and 
Prandtl's  criteria  are  all  included  in  our  separation  conditions. 

(3)  For  separated  flows  with  one  zero  u  line,  the  conditions 
set  for  FUS  and  zero  vortex  criteria  for  separation  are  ample  but 
not  necessary.  Thus,  these  flows  can  not  represent  all  kinds  of 
separated  flows.  For  separated  flows  with  two  zero  u  lines,  the 
conditions  set  for  KRS  and  zero  vortex  criteria  are  necessary  but 
not  ample.  The  points  where  these  conditions  are  satisfied  are  not 
necessarily  separation  points.  In  a  word,  KRS  and  zero  vortex 
criteria  are  not  universal.  When  calculating  steady  flow  over  a 
downstream-moving  wall,  many  workers  found  that  the  >RS  conditions 
held  in  most  cases.  This  is  because  flow  with  two  zero  u  lines 
appears  in  most  cases  where  MRS  conditions  are  necessary.  It  is  als 
found  that  when  steady  flow  over  an  up stream- moving  wall  was  cal- 


culated,  the  MRS  conditions  did  not  hold  in  most  cases  where  flows 
with  a  zero  u  line  most  likely  appeared  and  the  MIS  conditions  were 
not  necessary. 

(4)  In  the  case  of  separated  flow  described  in  the  boundary- 
layer  equation,  if  flow  shows  two  zero  u  lines  and  the  MRS  conditions 
hold,  the  separation  point  is  generally  a  Goldstein- type  singular 
point  (viscous  singular  point).  In  this  case,  MRS  and  singular  ooint 
separation  criteria  are  basically  equivalent.  However,  if  flow  shows 
one  zero  u  line,  and  the  MRS  conditions  hold,  inviscid  singular  ooint 
appears  with  the  slope  of  a  streamline  passing  through  the  separation 
point  being  infinity.  Therefore,  the  singular  point  separation  cri¬ 
terion  is  not  universal.  When  the  MRS  conditions  hold,  either  vis¬ 
cous  or  inviscid  singular  point  may  appear. 

(5)  O'Brien's  thinking  is  correct.  He  believes  that  when  the 
streamlines  passing  through  the  separation  point  are  not  normal  to 
each  other,  zero-vortex  conditions  are  not  satisfied. 

3.  THREE-DIMENSIONAL  STEADY  BLOWS 

Before  three-dimensional  steady  separated  flows  over  fixed 
walls  are  investigated,  separation  line  needs  to  be  defined.  In 
literature,  such  sketches  as  shown  in  Fig.  3  are  often  used  to  il¬ 
lustrate  the  separation  line  geometrically.  The  separation  line  is 
the  intersection  of  the  separated-flow  surface  and  the  material  sur¬ 
face,  and  fluid  at  both  sides  of  the  surface  flows  upward  along  the 
surface.  With  separation  line  defined  in  this  way,  it  is  not  cor¬ 
rect  to  use  zero  friction  stress  as  a  criterion.  Recently,  two 
prevalent  separation  criteria,  envelope  and  asymptote  limiting 
streamline  criteria  are  still  under  discussion.  The  necessary  con- 


ditions  proposed  by  Bai-Zelikovin  (44)  for  three-dimensional 
steady  boundary-layer  separation,  *»-gradp=o,  5,  •  gradp>  o  , 

where  tw  is  the  wall  friction  stress  vector;  p,  the  pressure;  and 
ue,  the  velocity  of  the  outside  boundary  layer,  have  not  yet  ver¬ 
ified  experimentally.  It  is  not  correct  either  that  the  line  where 
the  absolute  value  of  the  friction  stress  is  minimum  is  considered 
as  the  separation  line.  It  is,  therefore,  essential  to  further  in¬ 
vestigate  conditions  for  determining  the  three-dimensional  steady 
flow  separation  and  the  flow  patterns  near  the  separation  line. 

We  have  studied  this  subject  and  established  the  following  condition 
for  determining  flow  separation  and  the  separated  flow  attachment 
using  above  definition  of  the  separation  line: 


I 
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Fig.  S  Illustrating  the  separation  line 

1-  streamline;  2-material  surface;  3-  separated  flow  surface; 
4-separation  line;  5-limiting  streamline;  6-material  surface; 
7-separated  flow  surface;  8-separation  line;  9-limiting  stream 
line;  10-streamline. 


(1)  Assuming  that  xyz  is  the  coordinate  system  for  the  normal 
curves,  where  y  axis  is  the  separation  line;  and  z  axis,  the  normal 
of  the  material  surface,  when  flow  separation  appears,  the  following 
conditions  at  any  point  ”0"  on  the  separation  line  are  satisfied: 

(*»/**),-  0  (**«/a*d*)t<  0  (3.1) 

where  u  is  the  velocity  component  in  x  direction  (  the  normal  of 
the  separation  line).  Obviously,  equation  (3.1)  is  an  extension  of 
two-dimensional  Prandtl's  criterion. 

(2)  When  the  separated  flow  attaches  the  material  surface, 
the  following  conditions  hold  at  any  point  "o"  on  the  attachment 
line: 

0  <( 0  - •(*?/**)  (3.2) 

It  should  be  noted  that  conditions  (3.1)  and  (3.2)  hold  for  NS 
equation  and  for  the  flow  described  in  the  boundary-layer  equation. 

By  using  the  above  conditions,  we  can  draw  the  following  con¬ 
clusions: 

(1)  For  the  flows  described  in  NS  equation,  the  separation 
line  is  a  limiting  streamline;  The  limiting  streamlines  around  the 
separation  line  converge  upon  it  and  take  it  as  their  asymptote 
(Fig.  4(a)). 

(2)  For  the  flow  described  in  the  boundary-layer  equation, 
flow  generally  exhibits  a  Goldstein-type  singularity  in  the  normal 
direction  of  the  separation  line,  and  the  separation  line  is  the 
envelope  of  the  limiting  streamlines  (Fig.  4(c)). 

(3)  For  the  flow  described  in  NS  equation,  the  attachment 
line  is  a  limiting  streamline,  the  limiting  streamlines  in  its 
vicinity  taking  the  attachment  line  as  their  asymptote  and  generally 


Fig.  4  Limiting  streamline  patterns  near  separation  and  attachment  tines 
(  a  )  Separation  line  for  NS  eqnationi  (  b  )  Attachment  line  for 
NS  equation'  <  c  )  Separation  line  for  boundary  layer  equation. 


1 -limiting  streamline;  2-separation  line;  3-attachment  line, 
diverging  outward  (Pig.  4(b)). 

For  Newtonian  fluid,  since  the  limiting  streamlines  are  fric¬ 
tion  stress  lines,  the  above  drawn  conclusions  are  also  adapted  to  the 

friction  stress  lines. 

based  on  the  above  conclusions  and  by  solving  the  boundary 
layer  equation,  C.  K.  Wang,  et  al.  inevitably  arrived  at  the  con¬ 
clusion  that  the  separation  line  is  an  envelope  of  the  limiting 
streamlines.  3y  solving  NS  equation  as  the  starting  point,  however, 
it  can  be  concluded  that  the  separation  line  is  a  convergent  asymp¬ 
tote  of  the  limiting  streamlines  in  the  vicinity  of  the  separation 
line.  In  this  sense,  therefore,  Legendre  and  Lighthill's  ideas 
are  also  correct.  The  reason  for  this  is  that  the  boundary  layer 
approximation  results  in  a  singularity  of  the  solution  at  the  sepa¬ 
ration  line,  so  that  aberration  of  the  limiting  streamlines  occurs 


and  the  separation  line  which  should  originally  be  a  convergent 
asymptote  of  the  limiting  streamlines  under  a  real  condition  changes 


to  an  envelope  of  the  limiting  streamlines. 


In  addition,  we  have  discussed  the  origin  and  development  of 
the  separation  line  and  attachment  line,  and  pointed  out  that  the 
origin  of  the  separation  line  assumes  two  initial  states,  as  shown 
by  Mask ell.  and  K.  C.  Wang:  (l)  Being  originated  at  a  normal  point 
with  a  nonzero  friction  stress  component.  In  this  case,  the  sepa¬ 
ration  line  behind  the  origin  is  a  monotonic  smooth  curve  (Fig.  5 
(b)).  The  streamlines  ahead  can  pass  through  both  sides  of  the 
curve.  (2)  Being  originated  at  a  singular  point,  where  the  fric¬ 
tion  stress  component  is  zero.  In  this  case,  the  separation  line 
behind  the  origin  encircles  a  separation  region  (Fig.  5(a))  and 
separates  the  streamlines  ahead  from  the  region.  During  the  back¬ 
ward  development  of  the  separation  line,  this  line  extends  back¬ 
ward  until  a  singular  point  appears,  and  terminates  at  the  sin¬ 
gular  point  in  a  certain  state.  Based  on  the  singular  point  theory 
for  a  differential  equation,  it  can  be  verified  through  analyzing 
the  pattern  of  the  limiting  streamlines  in  the  vicinity  of  the  sin¬ 
gular  point  on  the  separation  line  that  if  the  origin  of  the  separation 
line  is  a  singular  point,  the  line  can  only  be  of  saddle-type,  and 
if  the  terminal  of  the  separation  line  is  a  singular  point,  the 
line  should  be  of  node-  or  focus- type. 

For  the  attachment  line,  the  corresponding  conclusions  can 
also  be  drawn:  if  the  origin  is  a  singular  point,  the  line  should 
be  of  node-  or  focus-type;  and  if  the  terminal  is  a  singular  point, 
the  line  should  be  of  saddle-type.  It  can  be  verified  that  focus 
on  the  separation  line  or  attachment  line  is  a  pressure  extreme 
point.  Obviously,  the  singular  point  is  closely  related  to  the 


separation  line  pattern.  It  is,  therefore,  significant  for  the 


separation  line  studies  to  investigate  the  distribution  of  singular 
points  on  a  body  surface.  At  the  beginning  of  1960’s,  Davey(46)  and 
Lighthill [22]  pointed  out  that  the  singular  point  distribution  of 
the  limiting  streamline  equation  obeys  the  following  topological 
rule:  on  the  surface  of  a  rr.onoconnection  body,  the  total  number  of 
the  nodes  and  focal  points  is  2  more  than  that  of  saddles.  Later, 
the  topological  rule  extended  from  flow  over  a  wall  to  that  beyond 
the  wall.  Smith(47),  Hunt(33),  Tobak,  Peake(36),  et  al.  have  point¬ 
ed  out  that  when  a  monoconnection  body  is  immersed  in  a  flow  which 
is  homogeneous  at  infinite  distance: 

(l)  Assuming  that  a  monoconnection  body  is  constructed  by 
three-dimensional  body,  B  and  a  flat  wall,  P,  and  that  P  is  either  an 
infinite  or  ring-like  surface,  then  on  surfaces  E  and  ?,  the  number 


of  singular  points  on  the  limiting  streamlines  satisfies 
C(ZAO-( 25)^,-  0 


(3.3) 


where  1-total  number  of  the  nodes  and  focal  points;  and  2-total 

number  of  saddles. 

On 

(2)„  the  two-dimensional  cross-section  cutting  a  three-dimen¬ 
sional  body,  the  number  of  the  singular  points  on  the  streamlines 


satisfies 


(  2iV  +  -|-  ZJV'  )  -  (  25+  — 1-25')=-  1 

2  —*■**«+« 

.  Hunts  ft  A.  A  A  Aft  «AA«  (ft*  T 
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(3.4) 


where  1-total  number  of  nodes  and  focal  points  beyond  the  outline  of 
the  material  surface;  2-total  number  of  nodes  and  focal  points  on 
the  outline  of  the  material  surface  (also  called  the  total  number 
of  seminodes  and  semifocuses) ;  3-total  number  of  saddles  beyond  the 


outline  of  the  material  surface;  and  4-total  number  of  the  saddles  on 
the  outline  of  the  material  surface  (also  called  total  number  of 
semisaddles) . 


where  the  meanings  of  N,  S,  N',  and  S'  are  the  same  as  mentioned 
above.  As  an  example.  Pig. 6  gives  the  distribution  of  the  fric¬ 
tion  force  lines  and  singular  points  of  the  flow  around  a  hemisphere- 
cylinder  erected  on  a  surface,  which  is  obtained  by  solving  NS 
equation(48) .  Pig.  6(a)  shows  the  distribution  of  the  friction 
force  lines  and  singular  points  on  a  surface.  In  this  case,  the  sin¬ 
gular  point  on  the  intersection  of  the  cylinder  and  the  surface  is 
a  semi-singular  point.  As  shown  in  the  figure,  2-Y=  4 .  2-V'=o.25=3, 
25'=  4,1  »  thus,  ^2iV  +  — 2iV'  25+  — 1~25'  j  =  ~l  f  which  satis¬ 

fies  topological  rule  (2).  Fig.  6(b)  is  the  distribution  of  the 
singular  points  on  the  hemisphere-cylinder  surface.  If  the 
hemisphere-cylinder,  B  and  the  bottom  surface,  P  are  considered  as 
an  integral  whole,  their  intersection  is  the  tie-line  of  the  inner 
points  on  the  surface.  As  shown  in  the  figure,  2-/V=9,25=9  , 

thus,  (2iV-2S)..,=  0 ,)  ,  which  satisfies  topological  rule  (l). 

4.  K.  G.  WANG'S  N.ETHOD  AND  RELATION  BETWEEN  TWO-DIMENSIONAL  UNSTEADY 
PLOW  AND  THREE- DI KEN  SION  AL  STEADY  FLOW 

According  to  K.  C.  Wang(l7),  if  the  time  variable,  t  for  two- 
dimensional  flow  is  treated  as  a  variable  of  space,  a  three-dimen-, 
sional  flow  can  be  obtained.  The  equation  describing  the  wall 
limiting  streamlines  for  such  a  steady  flow  is 

dxdt=d  u/az),.oA2  (4.1) 

where  z  is  the  coordinate  system  perpendicular  to  the  material  sur¬ 
face;  u,  the  velocity  component  on  the  wall  in  x  direction;  Az,  the 
step  of  the  calculation  mesh.  K.  C.  Wang  obtained  the  pattern  of 
the  wall  limiting  streamlines  using  equation  (4.1)  with  the  numer¬ 
ical  solution  of  the  boundary  layer  equation,  and  calculated  the 


are  introduced,  where  u^  is  characteristic  velocity  of  the  flow 
field  and  L  is  the  characteristic  length.  If  pressure  is  a  func¬ 
tion  of  x  only,  then  dp/#**  0  .  Thus,  the  boundary-layer  coupled 
equations  for  two-dimensional  incompressible  unsteady  flow  can  be 


Except  that  the  initial  condition  is  y=0  (t=0),  the  corresponding 


boundary-layer  conditions  are 


u  *  to  «  0  ,  v 


Z  -**00  u 


V 


-A*  (4.4) 


By  transforming  equation  (4.2),  the  limiting  streamline 
equation  on  the  wall  is 


jix  _  (  »z  )oA*  _  jL  f  au  ) 
dy  u.  u.  V  a*  /, 


(4.5) 


It  is  not  difficult  to  verify  that  equation  (4.5)  is  the  same  as 
equation  (4.1). 

The  following  discussion  is  based  on  an  assumption  that  A  z  is 
very  small  (approaches  zero)  and  y  is  limited,  which  is  equivalent 
to  the  condition  where  t  is  very  large.  Using  the  same  technique 
in  section  3  to  investigate  the  transformed  three-dimensional 
steady  flow,  the  separation  conditions  (on  the  separation  line)  are 

(-£),“- r"8’  C 5Ksl"°  ).<0  (4-6) 

where  x  is  the  coordinate  perpendicular  to  the  separation  line; 
and  d  is  the  angle  between  the  separation  line  and  x  axis. 

According  to  the  separation  conditions  (4.6),  we  will  easily 
obtain  the  pattern  of  the  streamlines  near  the  separation  line,  as 
shown  in  Fig.  7. 

On  the  other  hand,  for  two-dimensional  unsteady  flow,  sepa¬ 
ration  conditions  (1)  and  (2)  are  given  in  section  2.  When^lz  is 
very  small  (approaches  zero),  it  can  be  concluded  that  these  two 

conditions  are  just  equivalent  to  (4.6),  and  Jjjp-Ctga  = - L  hy 

can  be  obtained,  where  h  is  the  distance  from  the  separation  point 
to  the  material  surface.  This  shows  that  according  to  K.  C.  Wang, 
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the  above-given  separation  conditions  for  two-dimensional  unsteady 
flow  and  three-dimensional  steady  flow  can  be  transformed  by  each 
other. 


5.  SUMMARY 

This  paper  reviews  the  recent  studies  on  the  criteria  for 
separation  of  two-dimensional  unsteady  flow  and  three-dimensional 
steady  flow  and  the  characteristic  features  of  flow  in  the  vicinity 
of  the  separation  point  and  the  separation  line,  points  out  the 
existing  problems  and  introduces  results  of  the  authors  work  on  the 
topic  to  date,  which  may  be  helpful  in  solving  the  existing  problems. 

In  order  to  further  investigate  the  universality  of  these  results, 
it  is  necessary  to  perform  numerical  experiments  on  flows. 

However,  the  scope  of  the  present  paper  is  only  limited  to  the 
neighbourhood  of  the  separation  point  in  the  case  of  two-dimensional 
unsteady  flow  and  in  the  vicinity  of  the  wall  separation  line  for 
three-dimensional  steady  flow.  Since  separated  flows  usually  in¬ 
volve  larger  flow  fields,  it  is,  therefore,  essential  to  further 
investigate  the  space  structure  of  the  separated  flows.  Finally, 
since  separated  flows  and  vortex  motion  are  related  to  each  other, 
an  understanding  of  this  relationship  is  of  significance  to 
engineering  calculations. 
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ABSTRACT 

The  longitudinal  stability  of  elastic  vehicles  is  studied  with  the  effect 
of  aer oelasticity  considered.  The  transfer  functions  of  an  elastic  vehicle 
with  the  first,  second  and  third  elastic  modes  taken  into  account  are  de¬ 
rived  from  the  equations  of  its  longitudinal  disturbance  motion.  The  fo¬ 
rmulas  obtained  can  be  easily  used  on  computers  and  extended  to  cases 
with  higher  elastic  modes. 

The  effect  of  aeroelasticity  on  the  longitudinal  stability  of  a  vehicle 
is  investigated,  and  a  method  for  choosing  the  proper  position  of  the 
gyro  is  proposed.  Finally,  the  aeroelastic  effect  on  the  values  of  aero¬ 
dynamic  derivatives  is  estimated. 

The  conclusions  drawn  from  the  analysis  are  as  follows, 

(  1  )  The  aeroelastic  effect  on  the  stability  of  a  vehicle  is  mainly  due 
to  the  effect  of  the  first  elastic  mode,  while  the  effects  of  higher  ones  are 
relatively  of  less  importance. 

(2)  From  the  viewpoint  of  securing  stability,  the  gyro  should  be  lo¬ 
cated  behind  the  valley  of  the  irt  elastic  mode  if  the  generalized  force 
F.-*  and  control  moment  Nl,  for  unit  control  deflection  possess  the  same 
sign,  and  ahead  of  the  valley  if  F,»  and  Ms  possess  opposite  signs.  In  the 
case  of  the  investigated  vehicle  the  gyro  should  be  located  behind  the  va¬ 
lley  of  the  first  elastic  mode  and  ahead  of  the  valley  of  the  second  one. 

(  3  )  The  elastic  aerodynamic  coefficients  for  the  investigated  vehicle 
have  values  about  10—20%  different  from  those  of  the  corresponding  rigid 
vehicle. 
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NOMENCLATURE 


Zw  Aerodynamic  force  derivative  in  z  direction  caused  by  the 
normal  impact  speed. 

Zq  Aerodynamic  force  derivative  in  z  direction  caused  by  the 
angular  speed  of  pitch. 

Zfc  Aerodynamic  force  derivative  in  z  direction  caused  by 
deflection  of  the  control  mechanism. 

Z»ji  Aerodynamic  force  derivative  in  z  direction  caused  by  the 
generalized  displacement  of  the  ith  bending  mode. 

Z>jt-  Aerodynamic  force  derivative  in  z  direction  caused  by  the 
generalized  displacement  velocity  of  the  ith  bending  mode. 
Damping  coefficient  of  the  ith  bending  mode  structure. 

Fiq  Generalized  aerodynamic  force  derivative  of  the  ith  bending 
mode  in  z  direction  caused  by  the  angular  speed  of  pitch. 
Generalized  aerodynamic  force  derivative  of  the  ith  bending 
mode  in  z  direction  caused  by  the  normal  impact  speed. 
Generalized  aerodynamic  force  derivative  of  the  ith  bending 
mode  in  z  direction  caused  by  deflection  of  the  control 
mechanism. 

Fj^  Generalized  aerodynamic  force  derivative  of  the  ith  bending 
mode  in  z  rirection  caused  by  the  generalized  displacement 
of  the  ith  bending  mode. 

Fj^.  Generalized  aerodynamic  force  derivative  of  the  ith  bending 
mode  in  z  direction  caused  by  the  generalized  displacement 
velocity  of  the  ith  bending  mode. 

?',w  Longitudinal  moment  derivative  caused  by  the  normal  impact 
speed . 
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Longitudinal  moment  derivative  caused  by  the  angular  speed 
of  pitch. 

Longitudinal  moment  derivative  caused  by  the  rate  of  the 
impact  speed. 

Longitudinal  moment  derivative  caused  by  the  generalized 
displacement  of  the  ith  bending  mode. 

Longitudinal  moment  derivative  caused  by  the  generalized 
displacement  velocity  of  the  ith  bending  mode. 

Longitudinal  moment  derivative  caused  by  the  generalized 
displacement  acceleration  of  the  ith  bending  mode. 
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1 .  INTRODUCTION  322 

In  the  longitudinal  stability  studies,  some  investigators  * 
have  considered  an  elastic  vehicle  as  a  structure  where  a  rigid 
vehicle  is  in  parallel  with  an  elastic  mode.  In  such  a  structure, 
however,  the  aerodynamic  coupling  effect  between  the  rigid  vehicle 
and  the  elastic  mode  is  neglected.  Other  workers  considered  this 
coupling  effect,  but  the  mathematical  expressions  derived  are  com¬ 
plicated  so  that  it  is  not  feasible  to  conduct  practical  calcu¬ 
lations  . 

In  the  present  paper,  the  aerodynamic  coupling  effect  between 
the  rigid  vehicle  and  the  elastic  mode  is  considered,  and  a  standard 
matrix  form  of  the  equation  of  disturbance  motion  with  parameters  , 
such  as  and  M*g,  is  presented,  which  makes  It  feasible  to  use  the 

standard  compute  code  during  calculation.  Compared  with  those 
introduced  in  literature  T3-6],  this  technique  has  the  potential  to 
consider  higher  elastic  modes. 
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The  effect  of  gyro  location  on  elastic  vehicle  stability 
is  analyzed  together  with  the  coupling  effects  of  the  elastic  mode 
and  the  control  system.  It  is  believed  that  the  suitability  of  the 
gyro  location  is  related  to  the  same  sign  of  the  generalized  force, 
Fis  and  the  control  moment,  ivs .  Results  of  universal  significance 
are  given.  The  effect  of  aeroelasticity  on  aerodynamic  force 
derivatives  for  a  rigid  vehicle  is  also  discussed. 

2.  EQUATIONS  OF  ELASTIC  VEHICLE  LONGITUDINAL  DISTURBANCE  NOTION 

Now  we  are  going  to  establish  equations  of  the  short-period 
vehicle  longitudinal  disturbance  motion  with  aeroelasticity  taken 
into  account. 

The  aerodynamic  force  and  moment  are  related  to  the  local 
angle  of  attack,  a  ( x ,  y,  t)  or  the  local  imoact  soeed, 
ui(*,  y,  i)  =  u0a(*,  y,  t ),  ,  and  the  local  angular  soeed  of 

pitch,  o<  *•  >'*  O  . 

Let  us  assume  x,  y,  z  as  the  coordinates  in  the  average  body 
axis  whose  origin  is  the  instantaneous  center  of  mass,  and  express 

IB 

elastic  deformation  as  £(*.  y *  <  )•  2  ♦/(*>  O*  ,  where  ♦/(*»  y) 

1*1 

is  mode  function,  and  7Jj.(t)  is  the  generalized  coordinate. 

The  equations  of  the  nth  mode  short-period  small  nerturbation 
can  be  expressed  as 

*(  O-J,  O  dxdy 

+  LLl^«X'  y ’  0<ixdy+Zih(t) 

(2. la) 

Odxdy+  J,  y,  ,)  dxdy 

+  J\  j,  t)dxdy+MMt  ) 


(2.1b  ) 


‘i.(  *  '  )+®5V  * )  — ^T)  (2.1c 

where  Qj_(t)  is  the  generalized  force,  and  mj_(t)  is  the  generalized 
mass.  They  can  be  expressed  as 

L  3‘>Jxd3’  (2‘2a 

m.(  I  )  =  Hyfim(x'  y  * '  y)ydxdy  (?.2b 


The  integral  terms  and  the  generalized  force  term  in  eqns. 
(2.1)  are  functions  of  w(x,y,t)  and  t?(z,y,t).  When  aeroelasticity 
is  taken  into  account,  the  latter  two  functions  should  be 


w(x,  x,  <)“w(0+  Y  <M*'  yK-(<)-  Y  *•#(*»  >H(t) 

i  *  1  i*l 

(2, 3a) 

fl 

fl«fl(x,  X,  <)*«(<)-  Y  tU*’  X)t)<(0  (2.3  b  1 

i  *  1 

where  (’)  denotes  a  derivative  with  respect  to  t,  and  (’),  that 

«( * ,  x ,  t) 

with  respect  to  x.  Substituting  the  values  of  w(x,y,t)^and  re¬ 
arranging  them,  equations  of  the  elastic  vehicle  longitudinal  dis¬ 
turbance  are  obtained: 


w(  «  )-«ofl(  t)»Zm  •  W(  t )+  Z,e(  t)+  2  CZ„,n,(  OD+^»6(  O 

i  *  1 

(2.4a  ) 

n 

t(  1  t  t  )  +  t  )+  Y]  CA/„.»1,(  t  ) 

i  *  1 

+  A/ii%(  <  )  +  Af;*,ti,(  t  ))  +M,?>(  t  )  (2.4  b) 

*i(  /  )  +  2;.«>,n.(  t  )+*>’V,(  /  )  =  Fi*H'(  t  )+F„9(  t  ) 

H 

.  +  Y  CFi,/n/(  /  )+F,;(n,(  <))+F,»&(  O  (2.4  c) 

i-i 


where 


M *.*  M F ,.i  F lltt  F F i* F, » 


are  coefficients 


relating  to  the  elastic  modes: 


y)dxdy 

Zi‘-$J.Jw*dx-  y>dxdy-LL^x- y)dxdy 
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M'-LL{^*,lx-  y)~{x-i&+J£k)',^x’  y)\dxdy 


a*A/i 


ax  dy 


♦.(*,  y)dxdy 


f,«—  — ■  r  r  *)<*«** 

mi  J  y  J  *  ***> 


M 
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dxdy 


4.(*.  y)dxdy 


F'-r =3^  J  y  J  •  ^7+,<  *  •  '  »i(  *  •  *  ^ 
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<,2.5; 

(2.6) 

(2.7) 

(2.8) 
(2.9) 

(2.10) 

(2.11) 

(2.12) 

(2.13) 

(2.14) 


3.  CALCULATION  OF  ELASTIC  VEHICLE  TRANSFER  FUNCTIONS  TAKING  INTO 
ACCOUNT  THE  FIRST,  SECOND  AND  THIRD  NODES 

In  order  to  make  the  calculation  easier,  Laplacian  transform 
was  conducted  for  eqns  (2.4)  to  get  a  standard  matrix  form,  so  that 
the  matrix  compute  code  can  be  used  in  this  calculation. 

Let  us  assume  the  third  vibration  mode.  Then,  the  matrix  form 


of  the  coupled  linear  equations  is  obtained  as  follows: 


Z„—  i  u„  +  Z, 
Ml  Ml-  t 


7-  Z 

Af*;,  A/;, 

0  0 
(F1;1-2C,<o1-  *  )  (F,,,— «>!) 


(^*;.-2Ctffl*-  s  ) 


Ft,t-co» 
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where 


'  Af  » + M  ;,Z„ + M  •{  jF  lm  +  M  %tFlm + M 

'M ,+ A/i(n0  +  Z,)  +MilFl,+M'itFu+Mi}FM 

1«.1/;1  +  MiZ;I+Af-;1(F1;|-2t1a>1)+M;.Fa;1+A/;-.Fa;i 

-A/,1+3fiZ,1+.V?1(F„  -®J)+AfvtFt,  +3fy,Fa,/ 

.-A/;t+3/iZ;t+A/-i1FI;f+A/;1(F,;1-2U®t)+^r-iaFJiJ 

= A/,,  +  A/.Z„:  +  M\;  +  A/v^F,,,-  a>i)  +  A/;/,,, 

-A/;|+A/iZ;}+A/vIF,;J+MwtF1;1+3/iJ(Fava-2Ca<oa) 

-M,$+A/iZ,J+.V;IF„,+M,!lFl,1+MvJ(Fa,,-a>i) 

;  Mh  +  MjZ..  +  M  %  ,F  lt+ MntFu+ M*  tF  u 


When  higher  vibration  modes  are  considered,  the  same  method  can  be 
used.  From  eqn.  (3.1),  we  have  the  transfer  functions: 


,  .  -4i_*  C,lHc/+t,J,+c/_+C,J,+C,!l_+C,l  +  C,_ _ 

A  *  S  l. i‘  +  OjS7  +  a,**  +  oaJ*  +  o,J4  +  oas*  +  OjS1  +  a7s  +  o$) 


(3.3) 


(3.4) 


,  .  4m  dtts7  +  dus*  +  dttsi  +  dlts*  +  du5i~dttst  +  dlts  +  dn 

a  “  s  (i,  +  a1s7  +  <j;J*  +  <v,  +  a4s4^';'«iS  +  <VI  +  0ii~a«) 
A,.  dt„s7  +  dzlsf  +  d..s « 4-  di;s4  +  dus*  -  dtls2  +  d.,s  +  d  „ 
A  =  *  (s*  +  a1s7  +  j.s8  +  fl555  +  a45,^o!i1-:oli?+^s+^ 

A,s  dias7  +  duS*  +  d^  +  d^s4  f  -f-  d;ts*  +  d„s  +  W„ 

A  *  *  (s'  +  o^’  +  ats’  +  aj^  +  a.i'  +  ais’  +  ajj'  +  o^  +  a,) 


(3.5) 

(3.6) 


where  0,  ,  and  d„-Flt,  d.,-Fi>,  di{1^F^ 

In  the  case  of  a  closed-rina  system,  the  angle  to  which  the 
gyro  is  subjected  contains  not  only  $  ,  but  also  the  partial  deflec 
tion  angle  resulted  from  various  modes  at  the  point.  Thus,  the 
gyro-subjected  angle  ^  is 

*  +c-«Pl(  /  ):pi, +  (-<*>;(  I  ))’1.  +  C-<P',(  l  (3.7) 

where  ^  (  1 ),  <Pj(  / ),  <Pj(  / )  are  derivatives  of  various  vibraxion 
modes  of  the  gyro  at  point  1  with  respect  to  x.  Then,  transfer 
function,  (s)  is 

«%,( s  )=*f,-<p;(  [  )»%,(  s  )-<?',(  /  )»%,(  *  )-*;<"  /  )**%,(  * ) 

f.S7  + 1 , s*  +  its*  +  f  ,J4  +/,**  +  Ils'-  +  lts  +  l.l  (^.8) 

s  (s'  +  a ,s7  +  0,5*  +  a,s‘  +  a4$4  +  a,ss  +a,  j1  +  <^s  +  a.) 

where  the  coefficients  of  the  terms  in  the  numerator  are 

l,mCl- «PU  I  )</,!-' *',(  J  Ki-'Pi(  I  )d”  «  -  0  ~  7  (3.9) 

It  can  be  shown  that  in  a  closed-ring  system,  \V^t  and  W*  have  the 
same  nominator  but  different  numerators.  When  the  root  locus 
method  is  used,  as  we  will  discuss  shortly,  the  extreme  points  of 
Wtft  and  Vfi  are  the  same,  but  the  zero  point  distributions  are  dif¬ 
ferent.  Results  of  the  sample  calculation  selected  in  this  pacer 
are  shown  in  Table  1. 


Table  1 .  Roots  of  the  determinants 


4 

0 

-  3.546  ±34.777  i 

-  4.594  1 321 . 192  i 

-  31.471 ±  761 .227  i  i 

i 

-  20.6461  1177.87  i 

M 

-  13.51 

- 

-14.9091788  336  i" 

-  42  7901  749.657  i 

-  18.4991  1179.50  i 

0 

-  3.6751  34.354  i 

_ 

1 

-  31. 444  1  75'’. 876  i 

-  20.652  1  1178  08  i 

Snx 

0 

-  3. Ml  133.797  i 

|  -  3.327  1  317.153  i 

- 

-  20.708  1  1178.11  i 

Anj 

0 

-  1.887  ±  33.778  i 

i 

-  5.3901321  387  i 

-  31. 701  1  761. 038  i 

- 

4.  EFFECT  OF  ELASTIC  BODY  ON  THE  CLOSE-RING  SYS TEN  STABILITY 

The  stability  of  a  closed-ring  system  can  be  discussed  in  terms 
of  the  distribution  of  the  zero  points  and  extreme  points,  and  their 
changes  in  an  open-ring  system.  In  this  section,  only  changes  in  the 
zero  extreme  points  of  an  elastic  body  will  be  discussed.  It  is  as¬ 
sumed  that  the  gyro  and  the  automatic  machine  are  ideal  proportional 
links. 

As  mentioned  above,  having  considered  the  effect  of  the  par¬ 
tial  deflection,  W Slo  ,  which  is  given  by  eqn.  (3.8),  is  used  as  the 
elastic  body  transfer  function  in  the  closed-ring  system.  This 
function  has  exactly  the  same  nominator  as  in  ecn,  (3.3).  It  is 
independent  upon  the  location  of  the  gyro. 

The  coefficients  of  various  terms  of  the  numerator  in  eqn.  (3.8) 
are  determined  by  eqn.  (3.9)  where  the  coefficients  d„,  d„- 

are  influenced  by  aerodynamic  coupling.  It  is  more  important  that 
I  ).  *,(  1 ).  <*>J(  /T  are  closely  related  to  the  gyro  location  and 
may  have  opposite  signs  at  different  points.  Thus,  selection  of 
different  1  (the  gyro  location)  has  appreciable  influence  on  coef- 


ficients,  t0-ly. 

It  is  shown  from  the  sample  calculation  that  coefficients  Cj_ 
are  all  positive.  Although  d^,  d2i  and  d3^  can  assume  either 
positive  or  negative  values,  the  corresponding  coefficients  for  the 
same  vibration  mode  have  the  same  sign,  since  every  mode  itself  is 
stable. 

If  one  point  is  found  to  make  the  last  three  terms  on  the 
right  side  of  eqn.  (3.9) 

-<pj(  1  l  )</.,-< vi(  l  )rf„-  0 

then,  is  equal  to  W^.  This  means  that  the  elastic  vibration 

effect  is  completely  vanished.  It  is,  however,  very  difficult  to 
determine  this  point  accurately.  A  suitable  region  should  be  given 
here,  e.g.,  (_<?;(  j  )c/„-<p;(  I  )</.,— <Pj(  l  )dai’)>  0  .  For  instance,  a  point 

is  chosen  where  f*P»(  / )™ — 0.24*  fjC  O“0. 195,  <p5(  /  )*  — 0.46o  ,  then, 

-*»•!(_/  )d,i-<P't(  I  )d.j-  qpj(  /  )</„>  0  .  Also,  When  i*0~6;  ,  0 

(see  Table  2)  and  /7=c7=0.1460x  io11 

The  zero  points  determined  by  coefficients,  lj_  are 
-9.161,  —  1 .856  ±  136 .6«,  -19.17  ±  595 .8«,  —20.22  ±  1271* 

Compared  with  the  corresponding  zero  points  of  A?  in  Table  1, 
except  that  the  last  pair  of  zero  points  slightly  moved  left,  all 
the  others  moved  considerably  right.  The  last  pair  have  only 
slight  influence  on  stability  due  to  the  coupled  pair  formation 
with  corresponding  extreme  point. 

Since  cj_  has  the  same  sign  as  MJt  the  signs  of  d-jj_,  d2i  and 
d^i  are  the  same  as  those  of  F-jj,  ?2$>  and  F3S,  respectively.  In 
order  f or- C'P'JC  l  )dlL+^',(  1  +  /  )d^and  Cj_  to  have  the  same  sign,  the 

gyro  should  be  located  slightly  behind  the  valley  of  the  ith 


vibration  mode;  conversely,  it  should  be  located  slightly  ahead  of 
the  valley.  In  the  case  of  the  sample  calculation  described  in 
this  caper,  the  gyro  may  be  located  behind  the  valley  of  the  first 
mode,  but  ahead  of  that  of  the  second  one.  Since  the  third  rode 
has  slight  influence  (  jd^-jj  is  small),  it  was  not  taken  into  ac¬ 
count  . 


5.  ESTIMATION  OF  THE  EFFECT  OF  ELASTIC  VIBRATION  ON  AERODYNAMIC 
FORCE  AND  MOMENTS 

From  eqn.  (2.4)  and  the  aerodynamic  derivatives  of  a  rigid 
body,  the  following  relations  can  be  obtained  (where  subscript,  e 
denotes  elastic  body,  and  r,  rigid  body): 
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where  the  aerodynamic  derivatives  with  bars  (e.g.,  cj..  cj— )  stand 
for  the  means  of  the  aerodynamic  derivatives  in  various  integral 
intervals,  rescectively ,  i.e., 
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In  this  example  with  elastic  vibrations  taken  into  account,  the 
aerodynamic  derivatives  are  10-20%  different  from  those  of  rigid 
vehicles  (see  Jig.  l) 
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6.  CONCLUSIONS 

Based  on  the  above  discussion,  we  can  conclude  that 
(1)  When  the  effect  of  aeroelasticity  on  vehicle  motion  is 
considered,  we  can  draw  a  conclusion  which  agrees  with  the  char¬ 
acteristic  features  of  the  common  vehicles  in  small  and  medium 
sizes:  the  effect  of  the  first  elastic  mode  on  the  vehicle 
stability  is  significant  while  that  of  higher  ones  are  not. 


(2)  In  order  to  avoid  elastic  vibrations  lowering  the 

stability  of  the  overall  missile  system,  the  gyro  should  be  located 
behind  the  valley  of  the  ith  elastic  mode  if  Fit  and  have  the 

same  sign,  and  ahead  of  the  valley  if  Fjj  and  K$  assume  the  op¬ 
posite  signs.  This  is  a  conclusion  of  universal  significance.  In 
the  case  of  the  sample  calculation,  the  gyro  should  be  located  be¬ 
hind  the  valley  of  the  first  elastic  mode  and  ahead  of  the  valley 
of  the  second  one. 

(3)  Using  the  theoretical  expression  established  in  this 
work  for  the  ratio  of  the  aerodynamic  coefficient  of  an  elastic 
vehicle  to  that  of  a  rigid  one,  the  sample  calculation  shows  that 
the  elastic  aerodynamic  coefficients  are  10-20%  different  from  those 
of  the  corresponding  rigid  vehicle. 
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ABSTRACT 


A  three-dimensional  nonlinear  method  for  determining  the  aerodyna¬ 
mic  performance  of  a  jet  flap  wing  has  been  developed.  The  method  uti¬ 
lizes  a  system  of  bound  and  free  vortices  with  vortex  lattices  to  represent 
the  wing*  jetwake  and  tip  vortices.  The  kinematic  boundary  conditions  are 
satisfied  at  the  wing*  as  the  ordinary  vortex  lattice  theory  states.  But  the 
feature  of  the  present  method  lies  in  satisfaction  of  both  kinematic  and 
dynamic  conditions  and  consideration  of  vortex  interaction.  The  solutions 
obtained  include  chordwise  and  spanwise  load  distributions,  from  which 
sectional  and  total  aerodynamic  quantities  are  derived.  In  the  solution  the 
jet-wake  sheet  and  the  tip  vortex  shape  are  determined  simultaneously 
with  the  nonlinear  effects  taken  into  account.  The  solutions  obtained  by 
the  present  method  agree  well  with  other  analytical  results  and  available 
data  for  a  jet  wing  with  large  aspect  ratio.  The  present  method  has  no 
limitation  on  the  aspect  ratio  and  is  applicable  to  the  cases  with  large  jet 
blowing  angles  as  well.  The  partial  span  jet  blowing  and  the  disturbance 
of  main  wing  blowing  to  the  second  wing  are  also  treated  in  this  paper. 
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1.  INTRODUCTION 


It  is  a  well-known  principle  to  produce  strong  lift  by  using 

jet- flapped  wings  (JFWs).  Application  of  such  wings  pushes  forward 

their  aerodynamic  performance  studies.  In  two-dimensions,  Spence1 

has  established  the  JFW  shape  theory,  which  states  that  when  the 

jet  thickness  tends  towards  zero,  the  jet  momentum  is  still  kept  a 

limited  value  while  the  jet  mass  flow  rate  tends  towards  zero. 

Spence’s  hypothesis  has  been  accepted  by  most  investigators.  In 

cases  where  the  jet  structure  and  viscous  effect  need  not  to  be 

considered,  this  hypothesis  agrees  approximately  with  the  actual 

JFW’s  physical  properties,  and  usually  gives  good  results.  The  two 

representative  methods  which  can  be  taken  as  examples  are  Halsey's 

elementary  vortex  distribution  and  Sato's  conformed  map.  In  three- 

dimensional  cases,  various  methods  based  on  lifting-line  and 

lifting-surface  theories,  such  as  those  proposed  by  Mask ell  and 
4  *5 

Spence,  Las,  Shen  and  Lopez,  and  Wasson,  have  been  developed. 

The  elementary  vortex  distribution  method  given  in  literature  (6 J 
is  based  on  the  linearization  theory,  and  the  boundary  layer  con¬ 
ditions  for  the  wing  and  jet  are  satisfied  over  a  given  surface. 
Calculation  results  for  large  aspect-ratio  wings  using  this  method 
agree  well  with  experimental  data.  In  the  case  of  small  aspect- 
ratio  wings  or  large  jet  blowing  angles,  however,  the  nonlinear 
effects  should  be  taken  into  account. 

Big  progress  in  the  vortex  lattice  theory  for  determining 

7-Q 

aerodynamic  performance  of  wings  has  been  made  recently,  "  and 
several  methods  have  tfeen  developed  to  solve  problems  caused  by  a 
vortex  singularity.  Among  these  methods  i?  the  concept  of  so- 
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called  "limited-nuclear  vortex".  Nevertheless,  further  investi¬ 
gation  is  needed,  since  utilization  of  the  vortex  nuclear's  em¬ 
pirical  data  may  lead  to  a  decrease  of  calculation  accuracy. 

In  the  present  investigation,  the  vortex  lattice  method  was 
used  in  the  JFW's  study.  The  wing,  jet,  wake  and  tip  vortices 
were  represented  by  vortices.  The  kinematic  and  dynamic  con¬ 
ditions  were  satisfied  at  the  wing  and  the  jet  surface.  The  shares 
of  the  jet,  wake  and  tip  vortices  were  simultaneously  obtained  as 
part  of  the  solution  with  the  nonlinear  effects  taken  into 
account.  The  present  method  has  no  limitation  on  the  aspect  ratio 
and  the  jet  blowing  angle,  and  is  applicable  to  wings  with  any 
shapes  of  surface,  cross-section,  and  jet  momentum  distributions. 

2.  KETHOD  AND  BOUNDARY  CONDITIONS 

Let  us  use  the  hypothesis  of  an  incompressible  potential  flow 
without  taking  separation  and  viscous  effects  into  account,  and  use 
Laplacian  equation 

v*4>=  o  ( i ) 

where  <p  is  the  velocity  potential.  The  solution  is  expressed  as 
a  vortex  lattice  which  satisfies  this  equation. 

The  boundary  layer  conditions  required  by  equation  (l)  can 
be  divided  into  kinematic  and  dynamic  conditions.  The  kinematic 
condition  is  that  the  air  flow  is  tangential  to  the  material  sur¬ 
face.  If  the  jet  location  is  known,  the  air  flow  can  be  taken  as 
the  material  surface  and  the  kinematic  condition  has  been  satis¬ 
fied.  In  addition,  the  jet  should  satisfy  the  dynamic  condition 
that  the  pressure  difference,  to  which  the  jet  is  subjected,  p  is 


in  equilibrium  with  the  centrifugal  force 

&P-  -j/r  (2) 

where  j  is  the  jet  flow  rate  in  unit  span;  and  r  is  the  partial 
curvature  radius  of  the  jet.  Based  on  Spence's  hypothesis]  eon. 

(2)  can  be  written  as 

Yv  =  cCj2r 

where  y  is  the  vortex  intensity;  v,  mean  of  the  velocities  above 
and  beneath  the  jet  surface;  c,  the  chord  length  of  the  wing;  and 
Cp,  the  jet  momentum  coefficient. 

Let  us  use  the  Kutta's  conditions  to  make  the  solution 
unique.  In  the  absence  of  the  jet,  these  conditions  are  satisfied 
at  the  back  fringe.  In  the  presence  of  the  jet,  the  jet  at  far 
lower  reach  needs  to  be  parallel  to  the  freestream  flow  direction. 

In  addition,  it  is  assumed  that  there  is  no  transverse  momentum 
transfer  within  the  jet,  i.e.,  the  momentum  is  kept  constant  along 
a  jet  vortex  line.  This  hypothesis  appears  fairly  reasonable  when 
the  spanwise  variation  of  the  Cp  distribution  is  not  appreciable. 

It  is  shown  that  the  JFW  lifting  surface  theory  presented  in 
this  paper  is  different  from  other  common  theories  in  mixed  bounda¬ 
ry  layer  conditions  and  unknown  jet  location,  which  makes  the 
equation  solving  more  difficult. 

3.  SOLVING  PROCESS 

Vortices  were  used  to  represent  the  wing,  wake,  the  tin  vor¬ 
tices  and  the  jet  surface  (Fig.  l).  Here,  the  wing  was  devided  in¬ 
to  several  surface  elements,  each  being  expressed  by  a  hoof  vortex. 
The  bound  vortices  were  on  the  wing  and  the  free  vortices  constructed 
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the  wake  and  the  tip  vortices.  The 
spanwise  bound  vortices  were  located 
at  %  chord  length  of  the  correspond¬ 
ing  surface  elements  and  the  control 
points  at  3/4  chord  length.  Wherever 
a  jet  existed,  the  wake  was  neglect¬ 
ed.  The  jet  was  expressed  by  a  vor¬ 
tex  lattice  system,  the  control 
points  of  which  were  located  at  3/4 
chord  length.  The  thickness  effect 
was  not  considered  in  this  calcula¬ 
tion. 


Since  the  locations  of  the  jet,  wake  and  tip  vortices  were 
unknowns,  iteration  solving  methods  were  used.  Basically,  there 
are  two  iteration  methods:  (1)  Start  with  a  jet  shape  which  satis¬ 
fies  the  kinematic  condition,  and  then  find  the  jet  shape  which 
simultaneously  satisfies  both  dynamic  and  kinematic  conditions; 
and  (2)  Begin  with  a  jet  shape  which  satisfies  the  dynamic  condition, 
and  then  find  one  which  satisfies  simultaneously  both  conditions. 
Halsey  has  pointed  out  that  the  first  is  not  a  dependable  method 
and  it  is  difficult  to  get  a  convergence  from  this  method.  Our 
results  have  already  shown  this  disadvantage.  This  is  the  reason 
why  we  chose  the  second  method  in  this  study. 

During  calculation,  we  utilized  a  method  in  which  the  co- 

1  1 

efficients  were  to  be  determined  to  find  the  jet  share.  Results 
show  that  introduction  of  suitable  functions  can  sneed  ut>  the 


iteration  process.  It  is  believed  that  the  functions  to  be  se- 


lected  are  not  unique  and  functions  in  different  forms  can  be  used 
to  express  the  jet  shape  in  defferent  problems. 


4.  RESULTS  AND  DISCUSSION 

Calculations  were  performed  for  typical  JFWs  in  this  inves¬ 
tigation.  The  calculation  conditions  are  shown  in  Figs.  2-4. 

Figs.  2  and  3  give  the  spanwise  distribution  of  the  normal  force 
and  the  chordwise  distribution  of  pressure  at  half  wing  span. 

The  data  obtained  are  lO^o  higher  than  those  given  in  literature  [6i 
especially  for  the  front  and  back  fringes.  Further  study  is  needed 
to  find  the  reason  for  this  difference.  One  of  the  reason  is  oro- 
bably  neglect  of  considering  the  wake  at  the  location  where  there 
is  a  jet,  namely,  neglect  Of  the  controlling  effect  of  the  wake  on 
the  jet,  resulting  in  higher  calculation  values  of  the  normal  force 
Fig.  4  gives  Das'  experimental  results  which  are  in  between  the 
calculation  values  given  in  this  paper  and  those  in  literature  (6> 
Fig.  6  gives  results  on  partial  span  jet  blowing.  Neglecting  the 
wake  at  the  jet  and  its  vicinity,  a  simplified  inviscid  model 
results.  Fig.  6  shows  the  calculation  results  for  the  combination 
of  the  wing  and  the  tail,  indicating  that  the  influence  of  the  main 
wing  jet  on  the  tip  vortices  of  the  second  wing  (tail)  did  not 
fully  developed.  Results  for  the  case  of  a  large  jet-blowing  angle 
and  a  small  aspect  ratio  are  given  in  Figs.  7  and  8,  which  show  the 
calculation  capability  of  this  method. 
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Fig.  5  The  shape  of  jet  and  wake  on  partial  span  jet  blowing 
(a)  vertical  and  tide  view,  (6)  crota-section. 

1-partial  span  blowing. 


Fig  6  Interference  of  main  wing  jet  on  tail 
(0)  vertical  and  side  view(  (fc)  iross-section. 

1-main  wing;  and  2-tail. 
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Fig.  7  The  shape  of  jet  and  wing 
tip  vertex  at  fy-eo’dll 
(•>  vertical  and  side  view,  (6)  cross-section. 
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Fig  I  The  shape  of  jet  and  wing  tip 
vertex  for  a  small  aspect-ratio  wingtH) 
(«>  vertical  and  tide  view,  (ft)  cross-section. 
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NUMERICAL  COMPUTATION  OF 


TWOPKASE  TURBULENT  COMBUSTION  PROCESSES 


Zhou  Xiaoqing 

(Shengyang  Aeroengine  Institute) 


ABSTRACT 


The  mathematical  models  and  numerical  methods  are  presented  for 
computer  prediction  of  aeroengine  combustor  flows  and  characteristics. 
Among  these  models  the  droplet  turbulent  diffusion  model,  the  modified 
k~«  turbulence  model  and  the  hybrid  Arrhenius/EBU  model  are  described. 
The  Eulerian  scheme  is  used  to  construct  conservation  equations  of  both 
gas  phase  and  droplet  phase  flows.  The  self-compiled  computer  code 
GEMCHIP  programmed  by  author  is  provided,  which  incorporates  SIM¬ 
PLER  and  DROPLET  algorithms  for  gas  flow  and  droplet  flow  predic¬ 
tions,  respectively.  The  block  correction  technique  is  adopted  to  facilitate 
iteration  convergence.  The  predicted  results  for  three  typical  twophase 
flows  are  given  to  show  the  applicability  of  this  method. 
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1.  INTRODUCTION 


To  date,  research  and  development  on  aeroengine  combustors 
are  based  mainly  on  testing  while  making  semi-empirical  analysis 
subsidiary.  Great  advances  in  computation-combustion  science  have 
been  made  in  recent  years,  and  numerical  computation  has  become  an 
important  technique  for  combustor  research  and  designs. 

This  paper  presents  mathematical  models  for  twophase  tur¬ 
bulent  combustion  processes;  gives  the  conservation  equations  of 
various  governing  processes;  and  describes  the  algorithms  prooosed 

1  _  Tt, 

and  adopted  by  the  author,  which  are  embodied  in  a  numerical 
ellipse  compute  code,  GEKCHIP  for  multiphase  combustion,  heat 
exchange  and  inter-diffusion  flows.  This  compute  code  can  deal 
with  single-  and  two-phase  stable  flows  in  Cartesian  coordinates 
or  cylindrical  coordinates;  laminar  or  turbulent  flows;  cold  or 
chemical  reaction  flows;  and  pre-mixed  or  diffusion  flames. 

Results  of  sample  computation  under  flow  conditions  for  some 
combustors  are  given. 


2.  ?■  ATHEISTICAL  NCDELS  OF  THE  PHYSICAL  PROCESSES 
2.1.  Fog  model 

Rosin-Rammler  functions  are  often  used  to  calculate  the 


initial  distribution  of  drool et  size. 


,  d  N,  N,S  rat  +  4  )/SVU/,(  r,  V 

/fc'“  dr,  ~  r,m  TC(  t  +  1  r,m  ) 

v  <  (  r<  Y  f  rt(  t  +  0/53  V  81 

elp  (  V  r,m)  l  rc(  i  +  I  )/S 3  J  f 


where  Y.~x  is  density  of  tne  total  croolet  number;  S  and  t,  the 
fogging  parameters;  rq,  the  radius  of  droplets;  rqm,  volumetric 
mean  of  the  droplet  radius;  and  p,  Gama  function. 


In  the  numerical  computation,  n^^  is  used  to  replace  fn  l 


r  „  1 
'*  K*~ 
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2.2.  Droplet  Turbulent  Diffusion  Model 

Previous  computations  only  considered  the  orbital  motion  of 
droplets.  In  fact,  droplets  diffuse  turbulently  under  the  in¬ 
fluence  of  turbulent  air  mass.  The  following  droplet  turbulent 
diffusion  model  has  been  proposed  by  the  author,  based  on  Ward’s 


analytical  work. 


a-,=»K 
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where' <Tnx  and  (Tg  are  turbulent  Schmidt  number  of  droplet  and  that 
of  gas,  respectively;  k  and  6  are  the  turbulence  dynamic  energy 
and  dissipation  rate,  respectively;  and  Pi  and  Pg  are  the  density 
of  droplet  and  that  of  gas,  respectively. 

2.3.  Twophase  Turbulent  Model 

The  author  has  considered  the  turbulent  energy  exchange 
between  gas  phase  and  droplet  phase,  and  proposed  a  turbulent 
model,  based  on  k-e  relation.  In  the  model,  the  coupled  equations 


for  droplet  phase  turbulence  are 


—  V»i 
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where  6ni,  SVg  and  SVi  are  density  of  droplet  number,  the  tur¬ 
bulent  pulse  heights  of  the  gas  flow  rate  and  the  droplet  velocity, 
respectively.  The  ratio  of  the  droplet  turbulent  dynamic  energy  to 
that  of  the  gas  flow  is  k,/ k  ,  thus,  eqn.  k  can  be  revised 
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(7) 


where  <?*=tTjy^  . 

Eqn.  6  does  not  need  revision. 

2.4.  Droplet  Vaporization  and  Combustion  Kodel 

The  rate  of  the  droplet  size  change  during  vaporization  and 

combustion  is  ^rr= — — 

dt  Pi  c,  r, 

(8) 


where 


B  = 


N.=  1  +  0.276Re1'1•S'*/, 

<I>,  ==  ln(  1  +  B) 
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where  6-droplet  combustion;  and  7-droplet  vaporization. 

In  the  above  equation,  Yox  is  the  oxygen  mass  concentration;  Wox, 
the  molecular  weight;  vox,  the  chemical  equivalent  coefficient; 

I  and  q°  are  the  latent  heat  and  calorific  value  for  vaporization 
of  combustible  oil,  respectively;  and  Re  and  Sc  are  Reynolds  number 
and  Schm.it  numbex ,  respectively. 

2.5.  Gas  Phase  Turbulent  Combustion  Rod el 

Let  us  assume  that  combustion  is  a  single-step  chemical  reac¬ 
tion,  i.e., 


v,F  +  v0*0—  V.P 


(9) 


The  combustion  rate  is  the  smaller  one  between  the  rates  of 
the  gas  phase  reaction  process  (calculated  with  Arrhenius  equations) 
and  turbulent  mixing  process  (calculated  using  EBU  model) . 

Rr=  —  min( \Rr,itc  .  (  10) 
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where 


R',t*=-cM,k'  :[(^)  +(^t)  +(fe)  ] 

RftA,=  -BJWiYiYo*  exP(  -  E.  RJ,' 


where  E  is  the  activation  energy;  Ru  is  gas  constant;  and  Yp  is 
the  mass  concentration  of  the  combustible  material. 


3.  COUPLED  PARTIAL  DIFFERENTIAL  EQUATIONS  OF  THE  GOVERNING 
PROCESSES 


The  controllong  partial  differential  equations  for  two-phase 

turbulent  binary  chemical  reaction  flow  can  be  summed  up  as  the 

3 

following  two  forms. 


where  g  is  any  dimensionless  flow  variable;  Sj  and  S(JI/r  are  source 
of  action  terms  within  one  phase  and  between  two  phases;  subscript, 
i  is  the  coordinate  factor,  which  is  0  in  Cartesian  coordinates  and 
1  in  cylindrical  coordinates. 

Table  1  and  2  list  all  dimensionless  partial  differential 
equations  in  either  form  (11)  or  form  (12). 


4.  METHOD  OF  NUMERICAL  COMPUTATION 

Equations  ( 1 1 )  and  (12)  were  separated  into  finite  difference 
equations  using  controlling  volumetric  technique.  Interlocking 
grid  system  was  used  to  eliminate  the  fluctuation  of  the  pressure 
field.  Wall  function  technique  was  introduced  at  the  near-wall 
region  where  there  was  an  effect  of  molecular  viscosity.  As  shown 


Table. 2  Governing  equations  in  form  (12) 
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1-Type  of  equation;  2-Equation  of  droplet  momentum  in  x  direction; 

3-Equation  of  droplet  momentum  in  K  direction;  4-Eroplet  temperature  equation 
and  5- where. 


in  Fig.  1,  a  regular  grid  system  with  the  block-off  region  was 
used  to  simulate  the  irregular  contours  of  the  combustor  and  the 
stabilizer.  For  the  calculation  of  the  gas  phase  flow  field, 
SIMPLER  algorithm  was  used,  which  saves  1/3  to  1/2  time  as  com¬ 
pared  with  SIMPLE  algorithm,  and  relaxes  the  requirements  of  the 
subrelaxation  factor.*"  The  convection/diffusion  flux  on  the 
controlling  volumetrical  surface  was  computed  using  power  law 
program.  The  line-by-line  three-vertical-angle  matrix  algorithm 
capable  of  direction-changing  and  to-and-fro  scanning  was  used 
together  with  the  block  correction  technique  to  speed  up  the  con¬ 
vergence  and  improve  the  solution  stability. 

An  algorithm  called  LROPLET  was  developed  by  the  author  for 
droplet  flow  field  estimation.  This  algorithm  uses  power  law  for 
the  droplet  number  density,  and  wind- facing  method  to  calculated 
the  flux  at  the  intersection  plane  for  the  droplet  velocity  and 
temperature. 


Fig  1  Block -off  region  in  a  regular  grid  system 

1-nomial  boundary;  and  2-real  boundary. 


GEKCHIP  code  has  contained  the  above  algorithm.  The  code  i 
devided  into  two  part:  (1)  General;  and  (2)  For  the  User.  The 
first  part  lists  the  basic  algorithms  and  the  second  one  has 


corresponding  changes  based  on  the  type  of  problems,  geometrical 
shape,  boundary  and  entrance  conditions,  and  grid  arrangement. 

The  operation  sequence  of  GEKCHIP  is  as  follows:  firstly,  the  gas 
phase  flow  field  (without  droplets)  was  computed  using  SIMPLER, 
which  was  taken  as  the  primary  guessing  of  the  gas-chase  field 
for  two-phase  flow  calculation.  Secondly,  the  droplet  field  was 
computed  using  DROPLET.  Also  obtained  were  the  exchange  quan¬ 
tities  of  mass,  energy  and  momentum  between  gas  and  liquid  chases. 
These  construct  the  source  term,  S,Vr  in  the  gas-phase  equation 
and  can  be  substituted  back  to  SIMPLER  to  correct  gas  chase  field. 
Then,  DROPLET  was  re-entered.  This  kind  of  alternation  was  per¬ 
formed  until  the  convergent  criterion  for  both  gas  and  oil  phase 
flow  fields  was  satisfied. 

5.  APPLICATIONS 

GEKCHIP  code  has  been  applied  to  various  cases,  such  as 
single-phase  and  binary-phase  aerodynamics,  heat  transfer,  com¬ 
bustion  and  chemical  reactions.  In  this  section,  only  three 
samcles  will  be  given  due  to  limited  space. 

5. 1 . Two-phase  Cold  Flow° 

Fig.  2  shows  the  flow  parameter  distributions  for  the  gas- 
phase  and  droplet  phase  in  the  axial  direction.  It  is  shown  that 
the  smaller  the  droplet  size,  the  smaller  the  velocity  difference 
between  droplets  and  gas  flow.  T'ig.  3  shows  the  medium-size  drop¬ 
let  volumetric  fraction  curves  for  the  two-chase  flow  of  a  stream 
of  combustible  oil  gushing  from  the  center  of  a  truncated  cone. 
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Fig.  2  Twophuc  cold  flow  fields  ia  a  cylindrical  combustor 


(a)  Gas-phase  flow  rate  and  temperature  distribution  (dash  lines 
represent  single-phase  flow  and  solid  lines,  two-chase  flow) 

(b)  Droplet  flow  rate  distribution  in  the  axial  direction. 

(c)  Cylindrical  combustor.  1-air;  2-combustible  oil  and  air. 
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Fig.  3  Isograms  of  droplet  volumetric  fraction 
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Fig.  4  Tvopbasc  combustion  flow  fields  in  *  cylindrical  combustor 


(a)  Structure  of  the  flame  in  combustion  state  outside  the  dust 
fog. 

(b)  Concentration  and  temperature  distributions  in  the  axial 
direction. 

(c)  The  axial  distribution  of  the  combustion  efficiency  for 
various  nozzles. 

5.2.  Two-phase  Combustion  in  a  Cylindrical  Combustor^ 

Fig.  4(a)  shows  the  flame  structure.  The  flame  is  mainly 
located  outside  the  oil  cust  fog.  Fig.  4(b)  gives  the  axial 
distributions  of  the  combustible  material,  oxygen  concentration 
and  temperature  at  the  cross-section.  Fig.  4(c)  gives  the  axial 
distribution  of  combustion  efficiency  under  various  working  con- 
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Fig  5  Gu  pkut  hodograph  in  a  tvopbase  combustion  afterburner 

1-axis. 

5.3.  Two-phase  Combustion  in  the  Afterburner 

Fig.  5  gives  the  gas-phase  hodograph f or  the  afterburner  of 
a  small  culvert-ratio  fan  engine  with  a  parallel  inlet  gas  mixer/ 
pressurizer  and  a  stabilizer  with  single-row  v-shaped  slots.  The 
figure  clearly  shows  the  structure  of  the  return  air  region  end 
the  flow  acceleration  due  to  combustion.  It  is  shown  that  the 
influence  of  the  return  air  region  extended  to  the  lower  reach. 

6.  CONCLUSIONS 

The  method  and  GEKCHIP  code  presented  in  this  paper  have  been 
successfully  applied  to  several  typical  flows  in  combustors. 
Two-phase  ternary  code  is  expected  to  be  developed  from  this  work 
to  compute  performances  of  real  main  combustors  and  complicated 
afterburners.  A  close  combination  of  numerical  computation  with 
engineering  testing  and  semi-empirical  analysis  will  greatly  im¬ 
prove  research  and  development  on  combustors. 
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ABSTRACT 

An  analytical  method  of  gust  alleviation  and  other  active  control  pur¬ 
poses  (relaxed  static  stability  etc.  )  for  an  airplane  is  described.  Gust 
alleviation  uses  direct  measurements  of  atmospheric  turbulence  to  yield  a 
feed-forward  or  open-loop  control  law.  It  is  relatively  insensitive  to  chan¬ 
ges  in  flight  conditions.  Thus  a  combined  control  law  can  be  obtained 
which  contains  feed-back  and  feed-forward  components. 

The  stability  of  the  gust  alleviation  system  and  the  performances  of 
other  active  control  purpose  can  be  satisfied  by  obtaining  a  linear  opti¬ 
mal  control  defined  by  a  quadratic  performance  index  Thus  the  performa¬ 
nce  index  is  combined,  and  the  optimal  feed-back  control  law  is  strai¬ 
ghtforward  and  well-known. 

An  application  of  the  combined  control  law  to  a  flexible  airplane  is 
presented.  Calculations  have  been  performed  to  obtain  the  normal  accele¬ 
ration  responses  of  the  airplane  to  the  turbulence  and  the  state  responses 
for  the  other  active  control  purpose. 


Received  Inarch  1 


INTRODUCTION 

One  of  the  subjects  which  need  to  be  studied  for  some  airplane 
is  gust  alleviation.  Considering  that  the  application  of  active 
control  technology  (ACT)  to  airplanes  is  usually  malfunctional, 
to  combine  them  economically  and  effectively  is  a  subject  worth 
considering.  In  this  paper,  it  is  suggested  to  alleviate  airplane 
response  to  gust  by  using  an  open-loop  system  with  direct  gust  mea¬ 
surement.  The  requirements  for  other  active  control  functions  can 
be  met  with  a  feed-back  system.  The  combined  control  laws  for 
combined  gust  alleviation  and  other  active  control  purposes  are 
discussed. 

One  of  the  difficult  problems  which  one  faces  in  the  airplane 
gust  alleviation  system  designs  is  the  sensitivity,  i.e.,  big 
changes  in  the  system  parameters  with  flight  soeed  and  altitude. 

The  advantage  of  an  open-loop  system  with  direct  gust  measurement 
is  its  low  sensitivity.  Though  the  open-loop  system  has  its 
shortcomings,  an  improved  performance  can  be  obtained  by  combining 
an  open-loop  with  a  closed-loop  in  the  system. 

1.  PARTIAL  STATE  ALLEVIATION  FOR  DIRECT  GUST  NEASURENENT 

The  airplane  motion  equation  under  the  assumption  of  linear¬ 
ization  and  small  perturbation  is 

i(  t  )—Fx(  t  )  +  Gu(  t  )  +  /»,(  <  )  ) 

where  x(t)  is  the  state  vector;  u(t),  input  vector;  wg(t),  gust 

disturbance;  and  F,  G,  and  J  are  nxn,  nxm,  and  nxn  constant  matrices 
In  order  to  eliminate  the  effect  of  gust  on  eirolane  when  wt(t)  was 


measured;  '  the  following  equation  should  be  satisfied: 

G»(  t  O- 0  (?) 

However,  matrix  G  is  usually  a  square  matrix  and  its  inverse  can 
not  be  obtained.  Thus,  the  control  law,  u(t)  can  not  be  solved 
directly.  We  transformed  matrix  F  in  eqn.  (l)  to  a  diagonal  matrix 
and  assumed  the  transformed  matrix  as  T,  then  x-Tx  .  Thus,  ecn. 
(1)  could  be  written  as 

i-TFTmlx+TGu+TJu>t~A.x+TGa+TJwt  (3) 

where  A  is  a  diagonal  matrix  with  complex  numbers  in  it.  For 
convenience,  transformation  x=Mx  was  performed  to  eliminate  com¬ 


plex  numbers  in  A,  where 
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In  (4),  3-complex  characteristic  root;  and  4-real  characteristic 
root. 

Then,  eqn.  (1)  was  changed  to 


k^MTFT-'M-'x+MTGu+.UTJw, 


(5) 


MTFT-'M-'=a=  1 
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(1).  Ref.  4  discusses  the  methods  for  the  gust  measurement  and 


the  measures  to  reduce  the  measurement  errors.  Results  of 
flight  tests  under  various  flight  conditions  are  given. 


where  <Ti  and  CC j_  are  the  real  part  and  the  imaginary  part  of  the 
ith  complex  characteristic  root.  Eq.  (6)  has  r  complex  character¬ 
istic  root  and  n-r  real  ones.  Choose  m  state  variables  x-j(t)  in 
x(t),  which  play  a  major  role  in  gust  alleviation.  The  rest  n-r 
state  variables  are  X2(t)]  Eqn.  (5)  can  be  written  as 
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is  similar  to  a,  the  only  difference  between  them 
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o  aJ 

being  the  different  arrangement. 


g; 

LG.. 


anc 


are  similar  to  IVTC- 


and  MTJ,  respectively,  but  the  arrangement  is  different.  If  the 
following  equation  can  be  satisfied 

CGJ  {  u  > +  [/,!  =  0  (8) 

where  G-j  is  an  mxm  square  matrix  and  it  is  not  singular,  the 
control  law  is 

<«>--c<?.r,c/l3<wl>  (9) 

The  airplane  motion  equation  can  be  divided  into 


<*,>  -  Ca,D  {*,> 

ih  =  c  it)  +  ( c  cg,  r*  c/,3)  w 


(in) 


fx-j  is  independant  of  gust  and  fxg!  is  influenced  by  the  gust 
disturbance.  This  makes  the  {xij  -related  part  of  the  gust  response 
alleviated . 


2.  ALLEVIATION  OF  THE  OUTPUT  VECTORS 

It  is  of  significance  to  alleviate  airplane’s  measurable  out¬ 


puts  (such  as  acceleration  )  .  The  airplane  motion  and  output 
equations  are 


x(  t  )~Fx(  t  )+G,u1(  t  )  +  Jwt(  t  ) 

^i(  *  )==^iJC(  t  )+Dlul(  t  )-rEwt(  i  )  j  (  ^  ) 

To  alleviate  output  vector  y-|(t),  the  following  equation  should 
be  satisfied: 

y,(  X  )  =  C,x(  /  )  +  Dsu,(  1  )  +  Ewt(  0=0  (  1?J 

Generally,  E-j  is  not  a  square  matrix  and  has  no  inverse.  However, 
the  generalized  inverse  can  be  found  based  on  the  matrix  geometrical 
normal  number  and  eqn.  (12)  can  be  approximately  satisfied,  that 
is,  making  the  following  equation 

||CI*+D1B1  +  £w,i|,0l=(C,x+D;u1  +  j:u>,)TQi(C1x+D1u1  +  £'iiO=  v^Q,y,  ( 1  3) 

reach  a  minimum  value.  In  this  case,  Q-j  is  a  symmetrical  and  nor¬ 
mal  weighing  matrix  which  shows  the  alleviation  weighing  effect  for 
various  elements  in  y-](t).  Find  the  partial  derivative  of  eqn.  (13) 
with  respect  to  ui  and  make  it  equal  to  zero;  we  get 

2DlQlDlul  +  2D[QtClx+  2DIQ,Ew,=  0  (14) 

Solve  u-|f  and  the  control  law  is 

«1=-(DJQID,)*‘£>I01(CAx+ElBf)  (15) 

If  Xigix,  in  eqn.  (13)  is  replaced  by  y\Qlyl+u\Rlul  *  the  pro¬ 
cess  of  finding  the  generalized  inverse  has  already  considered  the 
requirements  for  the  controlled  energy.  R-)  is  a  symmetrical  and 
normal  weighing  matrix  which  controls  the  energy  weighing  quantity 

O 

at  the  control  surface.  The  corresponding  control  law  is 

( DJQA,  +  R,  )-Z?TO,(C1x+ Ewt )  (16) 

Under  the  gust  measurement  conditions,  eqns.(lc)  and  ( 1 6 ) 
contain  feed-forward  component  (D\QlD+Rl)'lDT1Q.Ewl  and  feed-back 
component  (D\Q,Dl  +  RXY'  *  DTQtC,*)  .  The  latter  component  can  pos¬ 


sibly  make  the  system  unstable.  A  stabilization  term  can  be  added 


to  the  feed-back  component  using  the  extreme  point  arrangement 
or  optimum  regulator  to  ensure  the  stability  of  the  overall  system, 
as  shown  in  the  following  sample  calculation. 

When  combining  gust  alleviation  with  other  active  control 
purposes,  assume  the  airplane  motion  and  outout  equations  as 

x=Fx+Gtu.  y^CtX+D^  (17) 

where  U2  is  the  control  output  satisfying  other  major  active  con¬ 
trol  purposes,  which  may  involve  an  operating  surface  different 
from  that  of  gust  alleviation.  Other  active  control  performance 
indices  can  be  written  as  second  order  forms,  such  as 


V  ■=  min4-  f  (ylQ.y.  +  ulR.aJdt 

ill  *  J  o 


(18) 


where  Qg  and  R2  are  symmetrical  and  normal  weighing  matrices. 
Combining  the  performance  index, equation  (18)  with  the  feed-back 
component  which  ensures  the  stability  of  the  gust  measurement  sys¬ 
tem,  the  combined  performance  index  is  written  as 


where 


V'  »minyJo  ( yTQy+uTRu)dt 
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(19') 


while  n  and  m  are  the  weighing  factors  measuring  various  weighing 
matrices  of  gust  alleviation  and  other  active  control  purposes. 
The  combined  performance  index,  eqn.  (19)  satisfies  eqn.  (18)  and 
at  the  same  time,  the  feed-back  control  makes  fust  alleviation 
stable  under  the  conditions  where  0'TQ:-vi+  “i^iO  ip  minimum. 


This  is  because  vector  y  in  eqn.  (19)  contains  output  vector  y-). 
Since  linear  optimal  regulator  control  is  independent  of  the 


system's  initial  condition  and  uncontrollable  disturbance,  the 
output  equation  in  eqn.  (11)  is 

y.-C.x+A*.  (20) 

The  optimal  control  to  assume  functional  V'  an  extreme  value  is 

u  =  -R-'&Px  (21) 

where  matrix  P  satisfies  the  following  algebraical  equation 


where  G 


G  ■ 


pTp+Pp-PGR-'&P+Q*  0  I 

Since  there  is  cross  term  when  eqn.  (19)  is 


expressed  in  state  variables,  transformation  is  needed  to  cancel 
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this  term;  Thus,  R,  F  and  Q  in  eqn.  (22)  are  given  as 
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Combining  eqns.  (16)— (21)  to  get  combined  control  law 


/r  ao, a  -i  \  r 

|  =  q  J+cA  cxyp^ 


D\QXE 


where  the  first  term  consists  of  two  parts:  feed-back  component 
which  alleviates  output  y^  in  the  gust  disturbance,  i.e.,  the 
feed-bajk  component  in  eqn.  ( 1 6 ) ;  and  the  control  quantity  which 
makes  V'  assume  a  minimum  value.  Because  V'  contains  y^  and  y^, 
it  satisfies  performance  index,  V  and  ensures  stability  of  gust 
alleviation.  The  second  term  is  the  open-loop  gust  alleviation 


term  resulted  from  gust  measurement. 


3.  SAMPLE  CALCULATION 

Let  us  take  the  variables  in  the  airplane  rigid-body  long¬ 
itudinal  motion  equation  and  body's  elasticity  equation  as  its 
state  variables.  Assume  that  the  control  surface  is  the  movable 


horizontal  tail  surface  and  the  wing  flap.  If  only  the  vertical 
component  of  the  gust  is  considered,  the  airplane  motion  equation 


r— i  .54 

1.0 

0 

0  1 
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Fi|  1  The  response  of  oatpat  >■ 
(normal  acceleration) 
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Fif-  2  The  response  of  output 
(other  active  control  purpose) 

1- witii  control 

2- no  control 


Assume  gust  as  u(l-costpt),  and  calculate  the  plane's  response  to 
gust.  Results  are  shown  in  Figs.  1  and  2,  which  compare  two  dif¬ 
ferent  conditions:  under  control  and  without  control.  In  calcu¬ 
lation,  it  was  assumed  u=20  m/sec  and  w=6 . 2852. 

4.  CONCLUSIONS 

A  gust  measurement  technique  for  gust  alleviation  are  sug¬ 
gested  for  airplanes  with  various  active  control  purposes.  The 
combined  control  laws  given  in  this  paper  can  meet  the  requirement 
for  gust  alleviation  and  other  active  control  purposes.  Response 
calculation  was  performed  for  an  airplane,  based  on  combined  con¬ 
trol  laws.  Calculation  results  and  data  comparison  show  obvious 
advantages  of  these  combined  laws. 
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CALCULATION  OF  SUPERSONIC  FLOWS 
AROUND  A  THREE-DIMENSIONAL  WING  AND  A  WAISTED  BODY 
WITH  CHARACTERISTIC  METHOD  IN  STREAK  SURFACE  COO'RDTNATES 


by 

Wang  Baoyu 

(Beijing  Institute  of  Aeronautics  and  Astronautics) 


ABSTRACT 


According  to  the  exact  characteristic  surface  equation  of  three  dime¬ 
nsions  and  the  physical  relations  along  bicharacteristics  in  stream  surface 
coordinates  of  Ref.  1 ,  the  pressure  distributions  of  circula-elliptic  delta 
wing  and  waisted  body  are  calculated  with  "real"  characteristic  method. 
The  results  are  compared  with  those  of  Refs.  2,3,4  and  6. 
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1.  INTRODUCTION 

A 

In  1960’s,  Bulter  first  performed  numerical  calculations 

for  inviscid  supersonic  flow  around  a  circular- elliptic  wing  at 

Kaod=3.5.  Ke  simplified  the  equations  by  using  high  supersonic 

characteristic,  and  obtained  surface  pressure  distribuions.  His 

results  agree  well  with  the  experimental  data.  In  1974,  Walkden 
2  3 

and  Caine  ’  calculated  the  same  wing  using  stream  surface  coor¬ 
dinates,  and  found  discrepancy  between  Suiter's  results  and 
theirs.  The  characteristic  equation  and  relations  of  physical 
property  variations  along  characteristics  (three-dimensional) 
on  which  Ref.  2  and  3  were  based  are  all  approximate.  Since  the 
term  of  the  partial  derivative  of  the  dependent  variable  with 
respect  to  jf3  was  moved  to  the  right  side  of  the  equation  in  Ref. 

5,  the  characteristic  equation  obtained  was  not  three-dimensional, 
but  equivalent  to  two-dimensional.  Moreover,  Ref.  2  and  3  used 
two  relations  of  the  physical  property  variations  along  character¬ 
istics  to  replace  two  relations  in  the  coupled  aerodynamic  equation 
and  solved  these  equations  using  difference  method,  not  the  char¬ 
acteristic  method.  Such  an  equation  solving  method  is  still  oven 
to  question.  In  Ref.  2,  the  calculation  results  with  different 
steps  were  interpolated  to  obtain  a  zero  mesh.  Ref.  3  mentioned 
tha;  calculation  could  not  be  continued  near  the  front  fringe,  tut 
Ref.  2  did  not  mention  this. 

The  author  obtained  an  accurate  three-dimensional  character¬ 
istic  equation  in  stream  surface  coordinates  and  an  accurate  re¬ 
lation  of  the  physical  property  variations  along  four  character- 

’Z  *7  p  p  a 

istics  over  two  sets  of  strerr  surfaces,  ^  =C^  and  5  =C  . 


In  the 


present  study,  this  "real"  characteristic  method  was  used  to  solve 
the  flow  fields  of  circular- elliptic  delta  wing  and  waisted  body. 

The  author  also  obtained  a  general  characteristic  relation  for 
three-dimensional  separated  flow  in  the  arbitrary  curve  and  stream 
surface  coordinates?  Considering  the  relation  given  in  Ref.  1  is 
lengthy,  it  is  not  given  here. 

2.  CALCULATION  KETHOD 

The  calculation  method  used  in  this  study  is  the  character- 

1  1 

istic  method  choosing  |  =x  in  a  stream  surface  coordinate  system. 

Ill 

If  the  aerodynamic  properties  at  the  surface  where  \  =x  =C  have 
been  caculated,  various  dependent  variables  at  point  A,  which  is 
on  surfaces  4  =c>»  i  =c  and  4  +^4  ,  need  to  be  calculated.  The 
calculation  procedures  for  space  points  and  boundary  points  are 
as  follows: 

(l)  Space  points.  Let  us  take  four  points  on  the  intersections 
of  4’=C\  |!=C1 2 3  and  (Fig.  1)  (i.e.,  points  6,  7,  S  and  9  at 

1- surface  4^ =C 

2- surface 

3- surface  ^  =c\ 


which  the  aerodynamic  parameters  can  be  solved  by  interpolating 
values  at  point  1  and  points  2,  3,  4  and  5)  to  construct  a  mesh 
cone  (characteristic  cone)  passing  through  A  (using  one-by-one 
approximation  technique).  Then,  take  points  6,  7  and  9  as  one 
group  and  points  7,  8  and  9  as  another  to  solve  the  aercdynar.ic 
parameters  at  A,  (P  )<i.  'fthn.  {*?)*i  and 

Parameters  at  point  A  are  the  averages  of  the  corresponding  values 
from  each  group.  Thus,  if  the  aerodynamic  parameters  at  calcu¬ 
lation  points  1,  2,  3,  4  and  5  on  surface  4I=jf'=cii  are 

known,  various  parameters  at  point  A  on  surfaces  '£!=CS  ,  4*»C* 

fc’+Afc1,  can  be  solved. 

(2)  Boundary  points.  The  calculation  method  for  the  boun¬ 
dary  points  are  similar  to  that  for  the  space  points.  The  rela- 
2  3 

tion  between  t^  and  t^  at  A  is  known.  For  point  5  and  9  on  sur¬ 
face  '4‘=Cl  (Fig.  1),  use  point  6  and  7,  and  point  7  and  8  to 

? 

determine  p  and  t1  at  point  A  and  then  to  arrange  the  corres¬ 
ponding  tv/o  values. 

Calculation  was  also  conducted  for  a  wais.ted  body  (axial  sym¬ 
metry,  two-dimensional).  Similar  calculation  method  was  used. 

There  were  no  space  points  6  and  8.  The  unknowns  to  be  determined 

2 

were  only  p  and  t^;  There  was  only  boundary  point  7,  and  p  was  the 
only  parameter  to  be  determined. 

In  this  method,  C.F.L.  conditions  must  be  satisfied  in  order 
to  make  a | 1  small  enough,  otherwise  unstability  would  appear  during 
calculation.  Such  an  unstability  did  occur  in  our  axial  symmetri¬ 
cal  flow  calculations.  There  was  no  unstability  found  in  the 
delta  wing  calculation,  since  selected  has  been  small  enough. 


3.  CALCULATION  RESULTS  AND  DISCUSSION 


Jl)  circular-  elliptic  delta  wing  (c*=0°).  The  geometrical 
shape  of  this  ''subsonic”  front  fringe  wing  can  be  found  in  Ref.  3. 
After  assuming  the  steps  as  A41— 0.003 985,  All  =  o. 0063754.  19635  > 

the  surfac e  pressure  distributions  on  the  tangential  surfaces 

A 

x  =0.417  and  0.683  were  calculated.  Results  are  shown  in  Pigs.  2 

p  '  0.15C 

and  3.  o.is 
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Fit.  2  Surface  pressure  of  win*  at  =  0  417  F>*  3  Surface  pressure  of  win*  at  *•  =  0  683 


1-Ref.  2  (zero  mesh);  2-Ref.  3; 
3-Ref.  4,  test;  4-Ref.  4,  theory; 
5- this  work. 


1-this  work;  2-Ref.  2  (zero 
mesh);  3-Ref.  3;  4-Ref.  4, 
theory;  5-Ref.  4,  test. 


Generally,  these  results  agree  fairly  well  with  those  given  in 
Ref.  2  and  arp  slightly  lower  than  those  in  Ref.  3°  The  character¬ 
istic  equation  and  the  relations  of  the  physical  property  variation 
along  the  bicharacteristics  on  which  the  present  work  was  based 
are  correct  and  the  calculation  method  is  a  "real"  characteristic 
one;  while  the  relations  on  which  Ref.  2  and  3  were  based  are  aoT- 
proximate,  and  the  methods  used  are  not  "real”  characteristic  ones 
(but  approximate  relations  of  the  physical  property  variations  a- 
long  the  characteristics  were  used  ).  The  difference  due  to  this 
is  conceivable.  The  "zero  mesh"  in  Ref.  2  contains  "real" 


characteristic  method  which  eliminates  the  difference  from  the 


approximation  algorithm.  Thus,  a  good  agreement  between  the  re¬ 
sults  of  this  work  and  those  given  in  Ref.  2  is  obtained.  The  cal 
culation  with  different  steps  in  Ref.  2  shows  that  the  smaller  the 
step,  the  lower  the  pressure  (see  Fig.  3  and  5  therein).  From 
this,  it  can  be  understood  why  the  author’s  results  are  slightly 
smaller  than  those  presented  in  Ref.  3. 

(2)  The  axial  symmetrical  flow  around  a  waisted  body  (of=o°) . 
Compared  with  Ref.  6,  this  work  gives  smaller  values  of  pressure 

distribution.  They  are  also  smaller  than  those  given  in  Ref.  7 

1  1 

(Fig.  4).  Near  x  =0.4  and  in  x  =0. 7-1.0  region,  difference  in 
pressure  distribution  is  appreciable  (but  still  very  small) . 

The  starting  equations  in  this  work  and  Ref.  6  are  correct.  In 
Ref.  6,  the  algorithm  is  not  a  ’’real"  characteristic  one  and  there 
was  no  "zero  mesh"  obtained. 


094 

Fit-  4  Surface  pressure  and  shape  of  waisted  body 

1-this  work;  2-Ref.  6;  3-Ref.  9;  4-test,  Ref. 7;  5-shape  of  waisted 
body. 

The  author  thanks  Zhang  Qinan,  who  did  a  lot  of  work  at  the 
beginning  of  the  computation,  and  Qian  Yiji,  who  provided  the  com¬ 
pute  code  for  the  Taylor-Faccoll  solution  to  flow  around  a  cone. 
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AESTRACT 


Based  on  the  thermodynamic  principles,  a  mixed  exhaust  turbofan 
cycle  is  split  into  an  inner  “regenerative*  subcycle  and  an  outer  simple 
subcycle  (Fig.  2)  with  the  borrowed  work  (Eq.  2)  and  the  borrowed  heat 
(Eq.  3)  refunded  to  the  inner  subcycle.  The  performance  of  the  engine  is 
shown  to  be  composed  of  subcycles' performances  through  respective  fraction 
ratios  <*  and  3  as  in  Eqs.  4 — 10.  Fig.  3  gives  a  thorough  view  of  the 
energy  and  efficiency  relation  in  the  subcycles  as  well  as  for  the  whole 
engine.  The  illustrated  example  in  Table  1  indicates  that,  although  combi¬ 
ned  with  an  efficient  inner  regenerative  subcycle,  a  low  compression  atd 
low  heating  outer  subcycle  gives  a  thermal  efficiency  inferior  to  the 
corresponding  simple  jet  engine.  The  overall  efficiency,  thrust  and  specific 
fuel  consumption  are  better  only  due  to  big.i  propulsive  efficiency.  The 
respective  functions  of  two  subcycles  in  the  er.  gine  are  thus  elucidated. 
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1.  INTRODUCTION 

Since  1960’s,  turbofan  engines  have  been  extensively 

developed  due  to  their  low  fuel  consumption.  However,  only  a 

1-5 

little  analysis  of  their  principles  has  been  conducted  so 
far,  most  of  which  are  in  terms  of  the  overall  engine  perform¬ 
ance.  In  the  present  paper,  the  engines  are  solit  up  into  inner 
and  outer  subcycles,  based  on  the  principles  of  self-induced 
compression  work  and  heating,  and  examples  are  given  to  show  the 
practical  thermodynamic  performances  of  each  subcycle.  The  sim¬ 
plified  solutions  of  the  pressure  ratios  for  three  compressors 
with  optimum  major  performances  are  obtained,  based  on  the  fixed 
specific  heat  of  the  overall  engine.  The  results  can  be  used  as 
a  reference  for  analysis  and  designs. 


2.  SEPARATED  EXHAUST  TURBOFAN  ENGINES 

Though  a  dual-functional  process  can  be  expressed  in  a 
temperature-entropy  diagram  (Fig.  l),  there  is  no  heat  source  in 


Fig.  1  Separate  Exhaust  Turbofan  (State  points  except  0,  9  refer  to  total  parameters, 

same  in  Fig  2) 
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the  outer,  and  a  cycle  does  not  occur,  based  on  the  second  law 
of  thermodynamics.  Thus,  analysis  can  not  be  done  on  it. 

However,  the  propulsive  efficiency  can  be  divided  into  the  ratios 
of  propulsive  work  to  effective  work  for  the  inner  stream  and 
outer  stream. 


T>»- 
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1  +B 


w, 
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i  1  4-  B  ' 

W,  +  w 
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where  dw  and  are  the  effective  work  part  ratio  of  the  inner 
stream  of  flow,  and  that  of  the  outer  stream  of  flow,  resoect- 
ively. 

It  is  shown  from  the  sample  calculation  in  Table  1  that  since 
there  exists  flow  resistance  due  to  the  turbo  and  fan  when  the 
effective  work  is  transferred  toward  the  outer  function,  the 
thermal  efficiency  of  turbofan  engines  is  lower  than  that  of  the 
corresponding  turbo  engines.  .  The  reason  for  an  improvement  of 
their  efficiency,  propulsion  and  fuel  consumption  is  the  improved 
propulsive  efficiency  due  to  separated  stream  of  the  low-soeed  jet 


3.  NIXED  TURBOFAN  ENGINES 

An  approximate  analysis  was  conducted  based  on  the  constant- 
overall-pressure  mixing.  The  same  analysis  can  also  be  used  for 
varying-overall-pressure  mixing.  Fig.  2  shows  the  thermodynamic 
process  which  was  split  up  into  inner  and  outer  subcycles.  The 
outer  compression  "borrowed  work"  from  the  inner 

w to.  -  B M ,  _ , ,/  /II.-  ( 1 + 9r )  (  2 ) 

and  the  outer  also  "borrowed  heat"  from  the  inner  (i.e.,  borrowed 


fuel  Bgr1  ) 


a +*!>*■«/ 


where  o<  and  ^  are  the  corresponding  part  ratios.  Accordingly, 
the  unit  propulsion  and  unit  fuel  consumption  are 
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It  can  be  seen  from  the  sample  calculation  in  Table  1  that  in 
terms  of  thermal  efficiency,  regeneration  improved  the  inner  sub¬ 
cycle;  however,  low  pressure  and.  low  temperature  made  the  outer 
subcycle  worse,  so  that  the  engine  performed  only  slightly  better 
than  the  corresponding  turbo  jet  engine.  It  is  the  high  propulsive 
efficiency  which  made  the  overall  efficiency,  propulsion  (ex¬ 
pressed  by  R/G^)  and  fuel  consumption  superior  to  those  of  a  tur¬ 
bo  jet  engine.  Since  the  heat  and  kinetic  energy  of  the  gas  flow 
in  the  outer  function  was  not  fully  used,  the  unit  propulsion  of 
the  turbofan  engine,  Rciw  is  lower  than  that  of  a  turbo  jet. 

After  slight  adjustment,  the  subcycle  method  can  be  used  to 
outer-functional  afterburning,  mixing  afterburning  and  dual¬ 
functional  afterburning  engines.  However,  when  designing  and 


analyzing  the  optimum  parameters  for  a  turbofan  engine,  a  simple 

1 

way  is  to  take  the  overall  engine* s  performance  into  account.  " 


4.  SIMPLIFIED  ANALYTICAL  METHOD 


In  this  study,  the  specific  heat  was  assumed  for  engine's 
overall  process.  Then,  the  performance  parameters  were  solved. 

The  condition  for  dual  functional,  constant-overall-pressure 
mixing  is 
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Assume  *  t  /=«*  *  and  B' =BA j/nj,  C=T\/T%'\\n\,  t ,=o,  * 

then  we  have  the  fan  pressure  ratio 

*  -.J±+C+B'C)x-Cx* 
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(11) 


V.'hen  -aJ— °  ,  fi'C***+(2C/T,)*— 04-C+5'C)/t,=  o 

From  that,  was  solved  for  the  compressor  pressure  ratio 
with  the  maximum  fan  pressure  ratio  (which  means  an  optimum 
afterburning  propulsion  and  lowest  fuel  consumption) : 


..  1  +*/l  +  B'(i+C+B'C)-z,  (1?) 

'  WVT, 

If  the  overall  pressure  loss  in  the  combustor  and  the  difference 

between  the  compressor  efficiency  and  fan  efficiency  were  not 

taken  into  account  (  (Tr=1=xr  and  3*=3),  then,  the  equation  could 

(4) 

be  simplified  into  A.M.  Lyul'ka’s  equation'  ; 

After  dial  functional  mixing 

Tt— U±bTji  .r:+ gr?  rj-rrcn; ~r- _  1)/T,*+Jer* 

1+5  1+b  r+j 

_  r?+(5+i/rit)r?  t\ 

1+5  nj(  i  +  b  ) x 

^  -  —  ✓  ^ 

P  Q 

Assume  .  —  I  1  -*  ~-l  ,  then  tne  temperature  decrease  of 

the  tail  jet  was  A7\*,~>)„(.P-Q*;(  i - L.) 


'  ntXtP-r>yf--b+B'(l-lhtmU)  +  i[C)x  +  1/tt  ,„t,C 

**- (1+5'  +  1/C)x 
^ _ _ . 

G* 


The  unit  propulsion  was 


:) 


From 


d  j-  •*  =  o  in  eqn.  (13),  Xt-  of  the  compressor  pressure  ratio 


dx  - ^  '  A 

(»*)*  when  the  propulsion  was  maximum  was 


x‘-2 GiX>-(£^-  P  +G.-G&  y  +  2Gt-^-x-GlGi-^-=0  (U) 

After  test  calculation,  x  had  only  one  real  root  from  1  to  3.727 
(*;=ioo) 

The  fuel  consumotion  rate  was 

h  -  1 

^  rs-TTCl  +  C"?  *  -D/15D  C,  3600 
1  +B  \H,  R*. 

K  S 

rs-rig-i/W-ayni)*  36ooc, 

R^_  U+B)IH. 

3600  /  c.a  I  K-Sx 


(15) 


(1 +  BKH 


/ 

y  2; 


_  '(P-Qx)(xx-<?1x  +  G’.)  - 7  1 

1®* - Zi — 7TZ - V JV2gc,J  • 


xx-G,x 


d{  ) 
minimum,  (re*)c 


V.x 


ui„  i  =  o  in  eqn.  (15),  when  the  fuel  consumption  rate  was 


K  f 


x  was 


/(/»-(?*  )(*»-(?, x+S,)7  /(P-Qx)(xx-Glx+Gt)  V, 

V  x'-G'X  yy  x‘-G,x 

q(x  -f  )[*4-2<?,x*  -(gl~^s  P  +Gi—G,GS)  x*  +  2G.  g  x  -(7:(73  g]  ( 16) 

+  2(x‘  — Cjx)* 

-  0 

After  test  calculation  with  equation  of  higher  degree,  x  only  had 


one  real  root  from  1  to  3.727. 


The  sample  performance  calculations  are  given  in  Fig.  4  and 
Table  2.  In  the  case  of  accurate  solutions  with  varying  specific 


Table  2  Variation  of  Performance  with  Compressor  Pressure  Ratio. 


/-  <C,=  0.24  k  «1.4>  |lfc»V*0-7 


nr*  -0.925 


=  0.848 
<Vt  *  0.91 


1-fixed  specific  heat;  2-simplified  solution;  3-accurate  solution 
with  varying  specific  heat. 


ffl4 

Pig  4  Variation  of  Performance  with  Compresior  Pressure  Ratio 

1- simplified  solution;  and  2-accurate  solution. 


heat  listed  herein  for  comparison,  the  mechanical  efficiency  _ 
7}m=0.99,  fuel  consumption  g  ,  outer  functional  overall  oressure 
coefficient  0^=0. 99  and  tail  Jet  efficiency  7^=0.9725  were  con¬ 
sidered  and  the  cooling  air  coefficients  were  introduced:  in 
combustion,  Y=0.91;  in  the  turbo  with  compressor,  y^=0.9456  and 
with  fan,  >^=0.97;  In  order  to  take  into  account  the  variations  of 
compressor  and  turbo  efficiencies  with  pressure  ratios,  the  multi¬ 
variant  efficiencies  equivalent  to  that  in  Table  1,  ()|^)k=0.848 
and  (7|p)»p=0*91  were  assumed.  This  means  that  the  ordered  ef¬ 
ficiency  was  nearly  unchangable.  The  big  circle  and  cross  in 
Fig.  4  and  the  underlined  data  in  Table  2  are  the  optimum  para¬ 
meters.  Generally,  the  initial  data  are  similar  to  those  in  Table  1 
It  is  shown  that  the  simplified  analysis  can  not  only  show 
the  tendency  of  variations  but  also  quickly  estimate  the  oressure 
ratios  of  the  two  compressors  with  optimum  propulsion  and  fan 
pressure  ratios. 
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